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ABSTRACT 

We present the discovery of the Kepler-20 planetary system, which we ini- 
tially identified through the detection of five distinct periodic transit signals 
in the Kepler light curve of the host star 2MASSJ19104752+4220194. From 
high-resolution spectroscopy of the star, we find a stellar effective temperature 
Tgff = 5455 ± 100 K, a metallicity of [Fe/H]= 0.01 ± 0.04, and a surface grav- 
ity of loggf = 4.4 ± 0.1. We combine these estimates with an estimate of the 
stellar density derived from the transit fight curves to deduce a stellar mass of 

= 0.912 ± 0.034 Mq and a stellar radius of = 0.9441°;°^° Rq. For three 
of the transit signals, we demonstrate that our results strongly disfavor the pos- 
sibility that these result from astrophysical false positives. We accompfish this 
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by first identifying the subset of stellar blends that reproduce the precise shape 
of the light curve and then using the constraints on the presence of additional 
stars from high-angular resolution imaging, photometric colors, and the absence 
of a secondary component in our spectroscopic observations. We conclude that 
the planetary scenario is more likely than that of an astrophysical false posi- 
tive by a factor of 2 x 10^ (Kepler-20b), 1 x 10^ (Kepler-20c), and 1.1 x 10^ 
(Kepler-20d), sufficient to validate these objects as planetary companions. For 
Kepler-20c and Kepler-20d, the blend scenario is independently disfavored by 
the achromaticity of the transit: From Spitzer data gathered at 4.5 fim, we in- 
fer a ratio of the planetary to stellar radii of 0.075 ± 0.015 (Kepler-20c) and 
0.065 ± 0.011 (Kepler-20d), consistent with each of the depths measured in the 
Kepler optical bandpass. We determine the orbital periods and physical radii of 
the three confirmed planets to be 3.70 d and 1.911q;2i -R® for Kepler-20b, 10.85 d 
and 3.07^1^;^? for Kepler-20c, and 77.61 d and 2.75l|]:^^ R(b for Kepler-20d. 
From multi-epoch radial velocities, we determine the masses of Kepler-20b and 
Kepler-20c to be 8.7 ± 2.2 M® and 16.1 ± 3.5 M®, respectively, and we place an 
upper limit on the mass of Kepler-20d of 20.1 M® (2 a). 

Subject headings: planetary systems — stars: individual (Kepler-20, KIC 6850504, 
2MASSJ19104752+4220194) — eclipses 



1. Introduction 



Systems with multiple exoplanets, and transiting exoplanets, each bolster confidence 
in the reality of the planetary interpretation of the signals and offer distinct constraints on 
models of planet formation. 



The first extrasolar planets were found around a pulsar (jWolszczan fc Fraillll992l ). and it 
was the multi-planetary nature — in parti cular the gravitational perturbations between the 
planets ( iRasio et al.lll992l : IWolszczanlll994j ) — which solidified this outlandish claim. Around 



Sun-like stars as well, the origin of radial velocity signals continued to be questioned by some, 
at the time multiple planets were found around ups Andromeda ( jButler et al.l 119991 ). The 
orbital configuration of planets relative to each other has sh ed light on a host of physical 
processes, from smooth radial migration i nto resonant orbits (ILee &: Pealell2002l ) to chaotic 
scattering into secular eccentricity cycles (jMalhotrall2008l : iFord et al.ll2005l ). N ow with ever 
grow ing statistics of ever-smaller Doppler-detected planets in multiple systems ( iMayor et al. 
20111 ). the formation and early history of planetary systems continues to come into sharper 
focus. 
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Concurrently, transiting exoplanets have paid burgeoning dividends, starting with the 
definiti ve proof that Doppler signals were truly due to gas-giant planets orbiting in close-in 



orbits (jCharbonneau et al.l l2000l : iHenry et al.l |2000[). Trans it lightcurves offer precise geo- 



metrical constraints on the orbit of the planet (jWinnI |2010| ). such that radial velocity and 
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the density of the planet and hence point to its composition 
Millen 2011i) . Transiting configur a tions also enable f ollow-up measure- 



Knutson et al.l 120071 : iTriaud et al.l |2010|) which inform on 



the mechanisms of planetary formation, evolution, and even weather. 

These two research streams, multipla nets and transiting planets, came toget her for the 



first time with the discovery of Kepler-9 (iHolman et al 



covery was enabled by data from the Kepler Mission (iBorucki et al 



2010 



Torres et a. 



2010 



2011). This dis- 



Koch et al.ll2010l ) 



which is uniquely suited for such detections as it offers near-continuous high-p recision photo- 



metric monitoring of target stars. Based on the first 4 months of Kepler data, iBorucki et al. 



(120111) announced t he detection of 170 stars each with 2 or more candidate transiting planets; 
Steffen et al.l ( l2010l ) discussed in detail 5 systems each possessing multiple candidate tran- 
siting planets. A comparative analysi s of the popula t ion o f cand idates with mul t iple pl anets 
and single planets was published by iLatham et al.l ( l201ll ) , and iLissauer et al.l (l2011al ) dis- 
cussed the architecture and dynamics of the ensemble of candidate multi-planet systems. 

The path to confirming the planetary nature of such Kepler candidates is arduous. At 
present, three stars (in addition to Kepler-9) hosting multiple transiting candidates have 



been presented in detail and the plane tary nature of each o 



l ished: These systems are Kepler- 10 (iBatalha et al. 



2011 



the candidates has been estab- 



(ILissauer et al.ll2011bl ). and Kepler-18 (ICochran et al. 201lf )~ Transiting planets are most 



ressin et al.l l2011al ). Kepler- 11 



profitable when their masses can be determined directly from observation, either through 
radial velocity (RV) monitoring of the host star or by transit timing variations (TTVs), as 
was done for Kepler-9bc, Kepler-lOb, Kepler- 1 lb cdef, and Kepler-lSbcd. When neither the 
RV or TTV signals is detected, statistical arguments can be employed to show that the plan- 
etary hypothesis is far more likely than alternate scenarios (namely blends of several stars 
containing an eclipsing component), and this was the means by which Kepler-9d, Kepler-lOc, 
and Kepler-llg were all validated. While such work proves the existence of a planet and 
determines its radius, the mass and hence composition remain unknown save for speculation 
from theoretical considerations. 

This paper presents the discovery of a new system, Kepler-20, with five candidate tran- 
siting planets. We vahdate three of these by statistical argument; we then proceed to use 
RV measurements to determine the masses of two of these, and we place an upper limit on 
the mass of the third. We do not validate in this paper the remaining two signals (and hence 
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remain only candidates, albeit very interesting ones, owing to their diminutive sizes), rather 



the v alidation of these two remaining signals is addressed in a separate effort (IFressin et al. 



20121 ). The paper is structured as follows: In §2, we present our extraction of the Kepler 



light curve (§2.1), our modeling of these data (and RVs) to estimate the orbital and physical 
parameters of the planets and star (§2.2), as well as limits on the motion of the photocen- 
troid during transit (§2.3) and a study of the long-term astrophysical variability of the star 
from the Kepler light curve (§2.4). In §3 we present follow-up observations that we use to 
argue for the planetary interpretation, including high-resolution imaging (§3.1) and Spitzer 
photometry (§3.2), and the spectroscopy we use to characterize the star and determine the 
radial velocity signal (§3.3). In §4, we present our statistical analysis that validates the 
planetary nature of the three largest candidate planets in the system. In §5 we consider 
the dynamics of the system, and in §6 we discuss the constraints on the composition and 
formation history of the three planets. 



1.1. Nomenclature 

Throu ghout the course of the Kepler Mission, a given star is known by many different 



names see 



Borucki et al.ll201ll . for an explanation of Kepler naming conventions), and we 
pause here to explain the relationship of these names in the current context. The star that 
is the subject of this paper is located at a =19^10™47^52, 8 =+42°20'19".4 (J2000). It 
was already known as 2MASSJ19104752-I-4220194, and in the Kepler input catalog it was 
designated KIC 6850504. After the identification of candidate transiting planets it became 
a Kepler Object of Interest (KOI) and was further dubbed K00070, and it appeared as such 



in the list of candidates published by iBorucki et al.l ( 1201 ll ) . Some authors have elected to 



denote KOIs using a different nomenclature, in which case K00070 would be identified as 
KOI- 70. After the confirmation of the planetary nature of three of these candidates it was 
given its final moniker Kepler-20. This paper describes that process of confirmation, but 
for simplicity we refer to the star as Kepler-20 throughout. The three confirmed exoplanets 
were initially assigned KOI designations representing the chronological order in which the 
transiting signals were identified, but to avoid confusion we will refer to them henceforth 
by their Kepler-20 designations in which they are ordered by increasing orbital period P; 
Kepler-20b (K00070.02, P = 3.70 d), Kepler-20c (K00070.01, P = 10.85 d), and Kepler-20d 
(K00070.03, P = 77.61 d). We will refer to the two remaining candidates as K00070.04 
and K00070.05, but note (as described below) that the period of K00070.04 (P = 6.10 d) 
is intermediate between those of Kepler-20b and Kepler-20c, and the period of K00070.05 
(P = 19.58 d) is intermediate between those of Kepler-20c and Kepler-20d. 



- 6 - 



2. Kepler Photometry and Analysis 
2.1. Light Curve Extraction 



Kepler observations of Kepler-20 commenced UT 2009 May 13 with Quarter 1 (Ql), 
and the Kepler data that we describe here extend through UT 2011 March 14 correspond- 
ing to the end of Quarter 8 (Q8), resulting in near-continuous monitoring over a span of 
22.4 mon ths. The Kepler bandpass s pans 423 to 897 nm for which the response is greater 
than 5% (IVan Cleve &: Caldwel]||2009f) . This wavelength domain is roughly equivalent to the 
l^ + i?-band JKoch et allboioh Thes e observations have been reduced and calibrated by the 
Kepler pipeline ( iJenkins et al.ll2010al ). The Kepler pipeline produces calibrated light curves 
referred to as Simple Aperture Photometry (SAP) data in the Kepler archive, and this is the 
data product we used as the initial input for our analysis to determine the system parameters 
(see below). The pipeline provides time series with times in Barycentric Julian Days (BJD), 
and flux in photo-electrons per observation . The data were initially gathered at long cadence 
(ICaldwell et al.ll2010l : iJenkins et al.ll2010bl ) consisting of an integration time per data point 
of 29.426 minutes. After the identiflcation of candid ate transiting planet s in the data from 



Ql, the target was also observed at short cadence (IGilliland et al. 



20101 ) corresponding to 



an integration time of 1 minute for Q2— Q6. We elected to use the long cadence version of 
the entire Ql— Q8 time series for computational efficiency. There are 29,595 measurements 
in the Ql— Q8 time series. The upper panel of Figure [1] shows the raw Kepler Ql— Q8 light 
curve of Kepler-20. The data are available electronically from the Multi Mission Archive at 
the Space Telescope Science Institute (MAST) Web sit^. 



2.2. Derivation of System Parameters 



The five candidate transi ting planets that are the subject of the paper were identified 
by the procedure described in lBorucki et al.l (120111 ). Four of them (Kepler-20b, Kepler-20c, 
Kepler-20d, K00070.04) are listed in that paper, and K00070.05 was detected subsequently. 

We first cleaned the Ql— Q8 long-cadence Kepler SAP photometry of Kepler-20 of 
instrumental and long-term astrophysical variability not related to the planetary transits by 
fitting and removing a second-order polynomial to each contiguous photometric segment. We 
defined a segment to be a series of long-cadence observations that does not have a gap larger 
than five measurements (spanning at least 2.5 hours). In this process, we gave no statistical 
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Fig. 1. — Kepler light curve of Kepler-20 at a cadence of 30 minutes. Upper panel: The 
normahzed raw SAP hght curves for Ql— Q8. The star is positioned on one of four different 
detectors, depending upon the particular quarter, which results in the most obvious offsets 
that occur roughly 4 times per year. The other discontinuities are due to effects such as 
spacecraft safe-mode events and loss of fine pointing. Lower panel: The SAP light curve 
after removing instrumental and long-term astrophysical variability via polynomial fitting 
(see §2.2p . In both panels transits of Kepler-20b are marked in green, 20c in red, 20d in blue, 
K00070.04 in cyan and K00070.05 in magenta. 



- 8 - 



weight to observations that fell within a transit. We then normalized the corrected light 
curve by its median, and this cleaned light curve is displayed in the lower panel of Figure [1] 

We then modeled simultaneously both this cleaned photometric time series and the 
radial- velocity measurements (described below in §3.3[ and listed in Table |5]) to estimate the 
orbital and physical parameters of the star and its candidate planets. The free parameters in 
the fit were the mean stellar density and the radial velocity instrumental zero point 7, and 
7 parameters for each of the 5 planet candidates i = {Kepler-20b, Kepler-20c, Kepler-20d, 
K00070.04, K00070.05}, namely the epoch of center of transit Tq^i, the orbital period Pi, 
the impact parameter 6j, the ratio of the planetary and stellar radii (i?p/i?^)i, the radial 
velocity semi-amplitude Ki, and the two quantities, (e cosa;)i and (e sina;)j, relating the 
eccentricity and the argument of pericenter Wj. The ratios of the semi-major axes to the 
stellar radius, (a/i?^)j, were calculated from and the orbital periods Pi assuming e = and 
that ^ sum of the planet masses. (We note that our observations do not constrain the 
eccentricity, but we include it so that our error estimates of the other parameters are inflated 
to account for this possibility. Similarly, we are not able to detect the radial- velocity signals 
Ki for Kepler-20d, K00070.04, or K00070.05, but by including these parameters, we include 
any inflation these may imply for the uncertainties on the mass estimates of Kepler-20b and 
Kepler-20c, and the upper limit on the mass of Kepler-20d.) 



We computed each transit shape using the analytic formulae of iMandel &: Agoll (120021 ). 
We adopted a fourth-order noii - linear limb-darkening law with coefficients fixed to those 
presented by IClaret &: BloemenI (1201 ll ) for the Kepler bandpass using the parameters Tgfj , 
\ogg, and [Fe/H] determined from spectroscopy ( §3.31) . Our approach implicitly assumes 
that all 5 transit signals are due to planets orbiting Kepler-20; the validation of the 3 largest 
planets, Kepler-20b, Kepler-20c, and Kepler-20d, is presented in §4. Using the validation 
approach presented in §4, we are not able to validate the remaining two candidat es K00070.04 



and K 00070.05. Instead this difficult problem is deferred to a subsequent study (IFressin et al. 



2OI2I ). We further assumed that the planets followed non-interacting Keplerian orbits, and 
that the eccentricity of each planetary orbit was constrained to be less than the value at 
which it would cross the orbit of another planet, e < 0.396 (Kepler-20b), 0.319 (Kepler- 
20c), 0.601 (Kepler-20d), 0.283 (K00070.04), and 0.325 (K00070.05). Finally, we included 
an additional error term on the radial velocities (beyond those appearing in Table E]) with a 
typical amplitude of 2 ms~^, to assure that we were not underestimating the uncertainties 
on the radial velocities (and hence the planetary masses). 

We included a prior on p^ as follows. We matched the spectroscopically-estimated 
Teff, log (7, and [ Fe/H] and the c orresponding error est imates (see §3.31) to the Yonsei-Yale 
evolution tracks (lYi et al.l I2OOII : iDemarque et al.l |2004[ ) with a Markov-Chain Monte Carlo 
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(MCMC) routine to produce posterior distributions of stellar mass and stellar radius i?* 
which in turn we used to generate the prior on used for the determination of the best-fit 
model. The posterior distributions of {a/R^)i were obtained by calculating {a/Ri,)i for each 
element of the Markov chain. 

We identified the best-fit model by minimizing the statistic using a Levenberg- 
Marquardt algorithm. We estimated the uncertainties via the construction of a co-variance 
matrix (these results were also used below in the estimate of the width of the Gibbs sample 
for our MCMC analysis). We then adopted the best-fit model (and its estimated uncertain- 
ties) as the seed for an MCMC analysis to determine the posterior distributions of all the 
model parameters. We used a Gibbs sampler to identify new jump values of test parameter 
chains by drawing from a normal distribution. The width of the normal distribution for each 
fitted parameter was initially determined by the error estimates from the best-fit model. We 
generated 500 elements in the chain and then stopped to exa mine the s uccess rate, and then 



were scaled the normal distributions using Equation 8 from iGregoryl (120111 ). We repeated 
this process until the success rate for each parameter fell between 22—28%. We then held 
the width of the normal distribution fixed. 

To handle the large correlation bet ween the model parameters, we adopted a hybrid 



MCMC algorithm based on Section 3 of iGregoryl (120111 ). The routine works by randomly 



using a Gibbs sampler or a buffer of previously computed chain points to generate proposals 
for a jump to a new locat ion in the parameter space in a manner similar to the DE-MC 



algorithm ( iTer Brakkll2006l ). The addition of the buffer allows for a calculation of vectorized 
jumps, which permits efficient sampling of highly correlated parameter space. Once the 
proposals drawn from the buffer reached an acceptance rate that fell between 22 — 28%, we 
held the buffer fixed. With the widths of the Gibbs sampler and buffer contents stabilized, 
we generated 4 chains, each with 1,000,000 elements. The generation of 4 separate chains 
permitted us to test for convergence via a Gelman- Rubin test. 

We adopted the median of the posterior distribution of each parameter as our estimate 
of its value, and we assigned the uncertainties by identifying the adjacent ranges of each 
parameter that encompassed 34% of the values above and below the median. We estimated 
the parameter distributions for the planetary masses and radii by combining the stellar mass 
and radius distributions from the evolution track fits with the model parameter distributions. 
The parameter estimates for the star are listed in Table [H and the parameter estimates for 
the confirmed planets Kepler-20b, Kepler-20c, and Kepler-20d are listed in Table [2l The 
light curves phased to the times of transit of each planet or candidate are shown in Figure [21 
In Figure [3] we show the radial velocities phased to the orbital periods of Kepler-20b and 
Kepler-20c. 
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We note that the values of a/R^ for K00070.04 and K00070.05 were misstated in Table 1 
of iFressin et al.l ((2012) • The correct values are 14.4^]^'2 and 31. Slug's? respectively. The values 
of the stellar parameters were also stated incorrectly in the same table, and should read 
Teff = 5455 ± 100 K, logg = 4.4 ± 0.1, [Fe/H] = +0.01 ± 0.04, ^;s ini < 2 km s -\ an d 
L/Lq = 0.71]']o29. These changes have no effect on the conclusions of iFressin et al.l ( 20121 ). 
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Table 1. Parameters for the Star Kepler-20. 



Parameter Value Notes 

Right Ascension (J2000) 19'^10™47=.52 

Declination {J2000) +42° 20' 19". 4 

Kepler Magnitude 12.498 

r Magnitude 12.423 

Spectroscopically Determined Parameters 



Effective temperature T^- 


ff (K) 


5455 ± 100 


A 


Spectroscopic gravity log 


:g (cgs) 


4.4 ±0.1 


A 


Metallicity [Fe/H] 




0.01 ± 0.04 


A 


Mt. Wilson S- value 




0.183 ± 0.005 


A 


logi?HK 




-4.93 ±0.05 


A 


Projected rotation -usini 


(km s"-*-) 


< 2 


A 


Mean radial velocity (km s ^ ) 


-21.87 ±0.96 


A 


Radial Velocity Instrumental Zero Point 7 (ms^^) 


_o 54+0.68 


B 



Derived stellar properties 

MassM*(MQ) 0.912 ±0.034 C 

Radius (Rq) 0.944l°;E;^g C 

Density (cgs) 1.51 ±0.38 C 

Luminosity L* (Lq) O.TIIqm C 

Age (Gyr) S.sttj^ C 

Distance (pc) 290 ±30 C 



Note. — A: Based on the spectroscopic analysis f i|3.3|l . 
B: Not a physical parameter, reported here for completeness. 
C: Based on a comparison of stellar evolutionary tracks to constraints from the 
spectroscopically-determined parameters ( i|3.3ll and the transit durations ( i|2.2l l. 
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Fig. 2. — Kepler light curves with an observational cadence of 30 minutes (black points) of 
Kepler-20, phased to each of the periods of the 5 candidate transiting planets (only data in the 
vicinity of each phased transit are shown). Kepler-20b, 20c, 20d, K00070.04 and K00070.05 
are shown from top to bottom. Blue points with error bars show these measurements binned 
in phase in increments of 30 minutes. The red curve shows the global best-fit model (see 
§2.2), which includes smoothing to match this 30 minute cadence. 
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Fig. 3. — Upper Panel: Radial velocities of Kepler-20 after correcting for the best-fit am- 
plitudes of Kepler-20b, Kepler-20d, K00070.04, and K00070.05, leaving the effect of only 
Kepler-20c and plotted as a function of its orbital phase of Kepler-20c. Individual measure- 
ments as shown as gray points and these values binned in increments of 0.1 phase units are 
shown in blue. The phase coverage is extended by 0.25 phase units on either side to show 
data continuity, but it should be noted the values in these gray regions are plotted twice. 
The red curve is the best-fit model for the radial velocity variation of the star after the sub- 
traction of the effect of Kepler-20b, Kepler-20d, K00070.04, and K00070.05. Lower Panel: 
Same as above, but showing the radial velocities (in gray, with binned points in blue) and 
model (in red) of Kepler-20 after correcting for effect of Kepler-20c, Kepler-20d, K00070.04, 
and K00070.05, leaving the effect of only Kepler-20b and plotted as a function of its orbital 
phase. 
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Table 2. Physical and orbital parameters for Kepler-20b, Kepler-20c, and Kepler-20d 



Parameter 


Kepler-20b 


Kepler-20c 


Kepler-20d 


Notes 


Orbital period P (days) 


3 69612iq+° '"^'""'-*3 

O.U»U±^±t)_g QQQQQg^ 


10.854092 ± 0.000013 


77 61184+0-00015 
' '-^ll^^-0. 00037 


A 


Midtransit time Tq (BJD) 




2454971.60758 ± 0.00046 


2454997.7271+0 00}^ 


A 


Scaled semi-major axis a/Ri, 


10 29+°-3'^ 


21-11?:? 


78.31^:^ ■ 


A 


Scaled planet radius Rp/Ri, 


n ni8c:t:+0. 00026 


0.02975 ± 0.00032 


0.0267010-00046 


A 


Impact parameter b 


goo+0.025 


594+0-018 


n 5RS+0-041 


A 


e cos {u)) 


-0 004+° °33 


-0 097+0 054 


-0 002+0-025 


A 


e sin{u}) 


-0 021+° °2l 


-0 011+0 031 

U.U±1_Q Q22 


-0 007+0-025 

U.UU(_Q Q22 


A 


Orbital inclination i (deg) 


86.50l°:36 


88.3910-16 


89 570+0-043 


A 


Orbital eccentricity e 


< 0.32 


< 0.40 


< 0.60 


A 


Orbital semi-amplitude K (ms~^) 


q 7^-1-0.76 
'^-1.09 


. on-l-1.03 

4.»J_o.9g 


1 2+1-0 


A 


Mass Mp {M^) 


8-7+22 


le.ilg 7 


< 20.1 (2 a) 


B 


Radius Rp (-R®) 


1 Q1+0:i2 

j^-yj^_o.2i 


Q n7+0-20 
•^■'^'-0.31 


r, 7C- + 0.17 

^- '0-0. 30 


B 


Density pp (gcm~^) 


6.5+2? 


2 ni+0-85 


< 4.07 


B 


Orbital semi-major axis a (AU) 


04537+0-00054 


0.0930 ±0.0011 


S453+0-0041 


B 


Equilibrium temperature Tcq (K) 


1014 


713 


369 


C 



Note. — A; Based on the joint modeling ( ^2.2| l of the light curve and radial velocities, with eccentricities constrained to 
avoid orbit crossing. 

B: Based on the results from A and estimates of Mi, and/or Ri, from Table [T] 

C; Calculated assuming a Bond albedo of 0.5 and isotropic re-radiation of absorbed flux from the entire planetary surface. 
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2.3. Limits on Motion of Photocentroid 



While the analysis above provides the parameter estimates of the five planet candidates 
under the assumption that each are planets orbiting Kepler-20, it does not address the 
concern that some or all of these five candidates result instead from an astrophysical false 
positive (i.e. a blend of several stars within the Kepler photometric aperture, containing an 
eclipsing component). In ^we use the BLENDER method to demonstrate that this possibility 
is extremely unlikely for Kepler-20b, Kepler-20c, and Kepler-20d, and it is this BLENDER work 
that is the basis for our claim that each of these three objects are planets. Another means to 



Batalha et al. 


2010; 


Torres et al. 


2011; 


Ballard et al. 



201 ll ) of Kepler-20 during times of transit, which we discuss below. Although we do not use 
the results presented below as part of the BLENDER work, we include a description of it here 
as it provides an independent argument against the hypothesis that the photometric signals 
result from an astrophysical false positive and not from planetary companions to Kepler-20. 

We use two methods to examine the Kepler pixel data to evaluate the location of the 
photocenter and thus to search for astrophysical false positives: (1) the direct measurement 
of the source location via difference images, the PRF centroid method, and (2) the inference 
of the source location from photocenter motion associated with the transits, the flux weighted 
centroid method. In principle both techniques are similarly accurate, but in practice the flux 
weighted centroid technique is more sensitive to noise for low signal-to-noise ratio (SNR) 
transits. We use both techniques because they are both subject to biases due to various 
systematics, but respond to those systematics in different ways. 



In our difference image analysis (ITorres et al.ll201l[ ). we evaluate the difference between 
average in-transit pixel images and average out-of-transit images. In the absence of pixel- 
level systematics, the pixels with the highest flux in the difference image will form a star 
image at the location of the transiting object, with ampl itude equal to the depth of the 
transit. A fit of the Kepler pixel response function (PRF, iBryson et al.ll2010l ) to both the 
difference and out-of-transit images provides the offset of the transit source from Kepler-20. 
We measure difference images separately in each quarter, and estimate the transit source 
location as the robust uncertainty-weighted average of the quarterly results. 

We measure photocenter motion by computing the flux-weighted centroid of all pixels 
downlinked for Kepler-20, generating a centroid time series for row and column. We flt the 
modeled transit to the whitened centroid time series transformed into sky coordinates. We 
perform a flt for each quarter, and infer the source loca tion by scaling the diff erence of these 
two centroids by the inverse of the flux as described in iJenkins et al.l (j2010d ). 
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Both the difference image and photocenter motion methods are vulnerable to various 
systematics, which may bias the result. The PRF fit to the diffe rence and out-of-transit pixel 



images is biased by PRF errors described in lBryson et al.l ((20101) ■ The photocenter technique 



is biased by stars not being completely captured by the available pixels. These types of biases 
will vary from quarter to quarter. Both methods are vulnerable to crowding, depending 
on which pixels are downlinked, which varies from quarter to quarter. We ameliorate these 
biases by taking the uncertainty-weighted robust average of the source locations over available 
quarters. Because the biases of these difference image and photocenter motion techniques 
differ, we take agreement of the multi-quarter averages as evidence of that we have faithfully 
measured the source location of the transit signal. 

Table E] provides the offsets of the transit signal source from Kepler-20 averaged over 
Ql— Q7 for all five planet candidates. The quarterly measurements and averages for the 
PRF centroid method are shown in Figure |H All the average offsets are within 2 sigma of 
Kepler-20. 



Table 3. Offsets between Photocenter of Transit Signal and Kepler-20 



Candidate 


PRF Centroid [arcsec] 


Significance'' 


Flux-wcightcd Centroid [arcsec] 


Significance'' 


Kepler-20b 


0.071 ± 0.25 


0.29 


0.41 ± 0.24 


1.72 


Kepler-20c 


0.021 ± 0.17 


0.12 


0.072 ± 0.22 


0.32 


Kepler-20d 


0.69 ± 0.74 


0.92 


3.07 ± 1.96 


1.59 


K00070.04 


0.24 ±0.51 


0.47 


2.14 ± 1.79 


1.20 


K00070.05 


0.73 ±0.45 


1.62 


1.76 ± 1.91 


0.93 



''offset / uncertainty 
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Fig. 4. — Quarterly and average source locations using the difference image (PRF centroid) 
method for Kepler-20b (top row left), Kepler-20c (top row center), Kepler-20d (top row 
right), K00070.04 (bottom row left), and K00070.05 (bottom row right). The light crosses 
show the individual quarter measurements and the heavy crosses show the uncertainty- 
weighted robust average. The length of the crosses show the 1 a uncertainty of each mea- 
surement in RA and Dec. The circles have radius 3 a and are centered at the average 
measured source location. The location of Kepler-20 is shown by the red asterisk and la- 
beled with its KIC number. In all the panels, the offset between the measured source location 
and the target is less than 2 a. 
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2.4. Stellar Activity and Rotation from the Kepler Light Curve 

While the polynomial- fitting approach in §2.2l is well-suited to cleaning the time series for 
the transit analysis, it annihilates any long-term variability, such as that due to rotationally 
modulated star spots. Subsequent to the analysis in §2.21 w e obtained a version of the Kepler 



photometry from Ql— Q8 using the new pipeline PDC-MAP (I Jenkins et al.ll201ll ). which more 
effectively removes non-astrophysical systematics in the photometry while leaving the stellar 
variability intact. We used this PDC-MAP corrected light curve to evaluate the rotational 
period and stellar activity of the star. We computed a Lomb-Scargle periodogram, and 
found the highest power peak at a period of 25 days, with a lobe on that peak at around 
26 days. This peak is also accompanied by significant power at periods between 24 and 
32 days. The distribution of periodicities and in particular the lobed, broad appearance of 
the peak with the highest power is strongly reminiscent of the activity behavior of the Sun, 
where differential rotation is responsible for a range of periods from approximately 25 days 
at the equator to 34 days at the poles. Indeed, the amplitude of spot-related variability 
on Kepler-20 is very similar to that of the active Sun, as measured during the 2001 season 
by the SOHO Virgo instru ment. Using the S OHO light curves treated to resemble Kepler 



photometry (as described in lBasri et al.ll201ll ). we compared the amplitude of variability of 



Kepler-20 and the active Sun; the two light curves (and the Lomb-Scargle periodogram of 
each) are shown in Figure We found that Kepler-20 has spot modulation roughly 30% 
higher in amplitude than that of our Sun. Our estimate of the rotation period (above) for 
Kepler-20 is consistent with both its spectroscopically estimated vsini < 2 kms~^ and an 
estimate of the rotation period, 31 d, based on its Ca H and K emission log-R^K (see §3.3p . 
Together, the period and variability indicate that Kepler-20 is very similar to (but perhaps 
somewhat more active than) our own star. 
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Fig. 5. — Upper panel: PDC-MAP corrected lightcurve for Kepler-20 from Ql to Q8 (black 
points), showing Kepler-20's intrinsic stellar variability after removal of instrumental effects. 
Orange points show the lightcurve for the 2001 active sun from the SQH O Virgo instrumen t 
g + r observations for comparison (lightcurve prepared as described in iBasri et al.ll2010l ). 
Lower panel: Lomb-Scargle periodogram for each of the two curves appearing in the upper 
panel. 
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3. FoUow-Up Observations 

3.1. High-Resolution Imaging 

In order to place limits upon the presence of stars near the target that could be the 
source of one or more of the transit signals, we gathered high-resolution imaging of Kepler-20 
with three separate facilities: We obtained near-infrared adaptive optics images with both the 
Palomar Hale 5m telescope and the Lick Shane 3m telescope, and we gathered optical speckle 
observations with the Wisconsin Indiana Yale NOAO (WIYN) 3.5m telescope. Ultimately 
we used only the Palomar observations in our BLENDER analysis (§4) as these were the most 
constraining, but we describe all three sets of observations here for completeness. 



3.1.1. Adaptive Optics Imaging 

We obtained near-infrared adaptive optics imaging of Kepler-20 on the night of UT 2009 
September 09 with th e Palomar Hale 5m telescope and the PHAR O near-infrared c amera 



(IHayward et al.l 120011 ) behind the Palomar adaptive optics system (ITroy et al.l 120001 ) . We 
used PHARO, a 1024 x 1024 HgCdTe infrared array, in 25.1 mas/pixel mode yielding a field 
of view of 25". We gathered our observations in the J (Aq = 1.25/xm) filter. We collected the 
data in a standard 5-point quincunx dither pattern of 5" steps interlaced with an off-source 
(60" East and West) sky dither pattern. The integration time per source was 4.2 seconds at 
J. We acquired a total of 15 frames for a total on-source integration time of 63 seconds. The 
adaptive optics system guided on the primary target itself; the full width at half maximum 
(FWHM) of the central core of the resulting point spread function was 0"07. 

We further obtained near-infrared adaptive optics imaging on the nigh t of UT 2011 June 



17 with the Lick Shane 3m telescope and the IRCAL near-infrared camera (ILloyd et al.ll2000l ) 
behind the natural guide star adaptive optics system. IRCAL is 256 x 256 pixel PICNIC 
array with a plate scale of 75.6 mas/pixel, yielding a total field of view of 19"6. We gathered 
our observations using the (Aq = 2.145/im) filter, and, as with the Palomar observations, 
we used a standard 5-point dither pattern. The integration time per frame was 120 seconds; 
we acquired 10 frames for a total on-source integration time of 1200 seconds. The adaptive 
optics system guided on the primary target itself; the FWHM of the central core of the 
resulting point spread function were 0.79". The final coadded images at J and Kg are shown 
in Figure O 

In addition to Kepler-20, we detected two additional sources. The first source is 3"8 
to the northeast of the target and is fainter by AJ ^ 4.5 magnitudes and AKs ~ 4.2 
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magnitudes. The star has an infrared color of J — K/^ = 0.19 ± 0.02 ma g which yields an 
expected Kepler magnitude of Kp = 16.1 ± 0.2 mag (iHowell et al.ll2011a( ). A much fainter 
source, at 11" to the southeast, was detected only in the Palomar J data and is AJ = 8.5 
magnitudes faii iter than the primar y target. The fainter source, based upon expected stellar 
Kp — J colors ( Howell et al.ll2011al ). has an expected Kepler magnitude of Kp = 21.0 ± 0.7 



mag. 

No additional sources were detected in the imaging. We determined the sensitivity 
limits of the imaging by calculating the noise in concentric rings radiating out from the 
centroid position of the primary target, starting at one FWHM from the target with each 
ring stepped one FWHM from the previous ring. The 3 a limits of the J-band and i^'-band 
imaging were approximately 20 mag and 16 mag, respectively (see Figure [7]). The respective 
J-band and i^-band imaging limits are approximately 8.5 and 5.5 mag fainter than the 
target, corresponding to contrasts in the Kepler bandpass of approximately 9 mag and 6.5 
mag. 



3.1.2. Speckle Imaging 



We obtained speckle observations of Kepler-20 at the WIYN 3.5m telescope on two 
different nights, UT 2010 June 18 and UT 2010 September 17. We gathered both se ts of 
observations with the new dual channel speckle camera described in iHorch et al.l ( 1201 ll ) . On 
both nights the data consisted of 3 sets of 1000 exposures each with an individual exposure 
time of 40ms, with images gathered simultaneously in two filters. The data collection, 
reduction, and image reconstr uction process are described in the aforementioned paper as 
well as in iHowell et al.l (l2011bl ). and the latter presents details of the 2010 observing season 
of observations with the dual channel speckle camera for the Kepler follow up program. 

On both occasions our speckle imaging did not detect a companion to Kepler-20 within 
an annulus of 0.05 — 1.8 arcsec from the target. The September observation yielded detection 
limits of 3.82 (in V) and 3.54 (in R) magnitudes fainter than Kepler-20. The June observation 
yielded limits of 3.14 and 4.92 fainter in V and R respectively. Therefore we rule out the 
presence of a second star down to 3.82 magnitudes fainter in V and 4.92 magnitudes fainter 
in R over an angular distance of 0.05—1.8 arcseconds from Kepler-20. 
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Fig. 6. — Palomar J (left) and Lick Kg (right) adaptive optics images of Kepler-20. The 
top row displays a 10" x 10" field of view for the Palomar J image and a 20" x 20" for the 
Lick Kg image. The bottom row displays zoomed images highlighting the area immediately 
around the target . 
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Fig. 7. — Left panel: The sensitivity limits of the Palomar J-band adaptive optics imaging 
are plotted as a function of radial distance from the primary target. The filled circles and 
solid line represent the measured J-band limits; each filled circle represents one step in 
FWHM. The dashed line represents the d erived corresponding limits in the Kepler bandpass 
based upon the expected Kp — J colors (IHowell et al.ll2011al ). Right panel: Same as above, 
but showing the sensitivity limits of the Lick Kg-hand adaptive optics imaging. The dashed 
line is based upon the expected Kp — Kg colors. 
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3.2. Photometry with the Spitzer Space Telescope 



An essential difference between true planetary transits and astrophysical false positives 
resulting from blends of stars is that the depth of a planetary transit is achromatic (ne- 
glecting the small effect of stellar limb-darkening), whereas false positives are not (except 
in the unlikely case that the effective temperatures of the contributing stars are extremely 
similar). By providing infrared time series spanning times of transit, the Warm Spitzer Mis- 
sion has assisted in the validation of many transiting p lanet system s, including Kepler-1 



( iFressin et al.l l2011a[). Kepler- 14 ( iBuchhave et a 



201lh. KeD ler-18 f ICochranet al.l l201lh 



Kepler-19 (IBallard et al.l l2011f l. and CoRoT-7 flFressin et al.l l2011bl l We describe below 
our observations and analysis of Warm Spitzer data spanning transits of Kepler-20c and 
Kepler-20d, which provide independent support of their planetary nature. 



3.2.1. Observations and Extraction of the Warm Spitzer Time Series 



We used the I RAC camera (IFazio et al.l |2004| ) on board the Spitzer Space Telescope 



( IWerner et al.ll2004j ) to observe Kepler-20 spanning one transit of Kepler-20c and two transits 
of Kepler-20d. We gathered our observations at 4.5 /im as part of program ID 60028. The 
visits lasted 8.5 hours for Kepler-20c and 16.5 hours for both visits of Kepler-20d. We used 
the full-frame mode (256 x 256 pixels) with an exposure time of 12 s per image, which yielded 
2451 and 4643 images per visit for Kepler-20c and Kepler-20d, respectively. 

The method we used to produce photometric time series from the images is described 
in lDesert et al.l ( l2009l ). It consists of finding the centroid position of the stellar point spread 
function (PSF) and performing aperture photometry using a circular aperture on individ- 
ual Basic Calibrated Data (BCD) images delivered by the Spitzer archive. These files are 
corrected for dark current, fiat-fielding, detector non-linearity and converted into flux units. 
We converted the pixel intensities to electrons using the information given in the detector 
gain and exposure time provided in the image headers; this facilitates the evaluation of the 
photometric errors. We adjusted the size of the photometric aperture to yield the smallest 
errors; for these data the optimal aperture was found to have a radius of 3.0 pixels. We found 
that the transit depths and errors varied only weakly with the aperture radius for each of 
the light curves. We used a sliding median filter to identify and trim outliers that differed in 
flux or positions by greater than 5 a. We also discarded the first half-hour of observations, 
which are affected by a significant jitter before the telescope stabilizes. We estimated the 
background by fitting a Gaussian to the central region of the histogram of counts from the 
full array. Telescope pointing drift resulted in fluctuations of the stellar centroid position, 
which, in combination with intra-pixel sensitivity variations, produces systematic noise in 
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the raw light curves. A description o f this effect, know n as the pixel-phase effect, is given 
in the Spitze r /IRAC data handbook (Reach et al. I2OO6) an d is well known in exoplanetary 



studies (e.g. ICharbonneau et al.ll2005l : iKnutson et al.ll2008l ). After correction for this effect 



(see below) we found that the point-to-point scatter in the light curve indicated an achieved 
SNR of 220 per image, corresponding to 85% of the theoretical limit. 



3.2.2. Analysis of the Warm Spitzer Light Curves 



We modeled the time series using a model that was a product of two functions, one 
describing the transit shape and the other describing the variation of the d etector sensitivity 
with time and sub-pixel position, as described in iDesert et al.l teOllaT). Fo r the transit 
light curve model, we used the transit routine OCCULTNL from lMandel fc Agoll ( 120021 ). This 
function depends on the parameters {Rp/R^,)i, {a/Ri,)i, bi, and To^j, where i = {Kepler-20c, 
Kepler-20d}, the two candidate planets for which we gathered observations. The contribution 
of stellar limb-darkening at 4.5 /im is negligible given the low precision of our Warm Spitzer 
data and so we neglect this effect. We allow only {Rp/R^)i to vary in our analysis; the other 
parameters are set to the values derived from the analysis of the Kepler light curve (see 
Table [2]). Because of the possibility of transit-timing variations (see §5), we set the values of 
To^j to the values measured from Kepler for the particular event. Our model for the variation 
of the instrument response consists of a sum of a linear function of time and a quadratic 
function (with four parameters) of the x and y sub-pixel image position. We simultaneously 
fit the instrumental function and the transit shape for each individual visit. The errors on 
each photometric point were assumed to be identical, and were set to the root-mean-squared 
residuals to the initial best fit obtained. 

To obtain an estimate of the correlated and systematic error s in our mea s urena ents, 
we use the residual permutation bootstrap method as described in iDesert et al.l ( 120091 ). In 
this method, the residuals of the initial fit are shifted systematically and sequentially by 
one frame and added to the transit light curve model, which is then fit once again and 
the process is repeated. We assign the error bars to be the region containing 34% of the 
results above and 34% of the results below the median of the distributions, as described in 



Desert et al.l (l2011bl ). As we observed two transits of Kepler-20d we further evaluated the 
weighted mean of the transit depth for this candidate. In Table HJ we provide a summary 
of the Spitzer observations and report our estimates of the transit depths and uncertainties. 
In Figure El we plot both the raw and corrected time series for each candidate, and overplot 
the theoretical curve expected using the parameters estimated from the Kepler photometry 
(see below). 



- 26 - 



The adaptive optics images described in §3. 1.11 reveal the presence of a star adjacent to 
Kepler-20. This adjacent star is 4.5 magnitudes fainter in J-band that Kepler-20, and located 
at an angular separation of 3.8", which corresponds to 3.1 IRAC pixels. We tested whether 
the measured transit depths have to be corrected to take into account the contribution from 
this stellar companion. We computed the theoretical dilution factor by extrapolating the 
J-band measurements to the Spitzer bandpass at 4.5 fim. We estimate that 1.6% of the 
photons recorded during the observation come from the companion star. For a blend of two 
sources, the polluted transit depth would be d/{l + e), where d is the unblended transit 
depth, and e = 1.6%. Since the effect on d is well below our detection threshold, we conclude 
that the presence of the contaminant star near Kepler-20 does not affect our estimates of 
the transit depths. 

We calculate the transit shapes that would be expected from the parameters estimated 
from the Kepler photometry (Table [2]) and overplot these on the Spitzer time series in Fig- 
ure [H The depths we measure with Spitzer are in agreement with the depths expected from 
the Kepler-derived parameters at the 1 a level. Our Spitzer observations demonstrate that 
the transit signals of Kepler-20c and Kepler-20d are achromatic, as expected for planetary 
companions and in conflict with the expectation for most (but not all) astrophysical false 
positives resulting from blends of stars within the photometric aperture of Kepler. 



Table 4. Transit Depths at 4.5 /jm from Warm Spitzer 

Candidate AOR Name Date of Observation [UT] Data Number Time of Transit Center [BJD] Transit Depth(%) 



Kepler-20c r41165824 2010-12-05 2291 2455536.0209 

Kepler-20d r39437568 2010-09-24 4451 2455463.4022 

Kepler-20d r41164544 2010-12-10 4383 2455540.9925 

Kepler-20d weighted mean — — — 



0.075 ± 0.015 

0.067 ± 0.016 
0.065 ± 0.011 
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Fig. 8. — Warm Spitzer transit light curves of Kepler-20 observed in the IRAC band-pass at 
4.5 fim spanning times of transit of Kepler-20c [upper half of figure) and Kepler-20d [lower 
half of figure). For each candidate, the raw and unbinned time series are shown in the upper 
panels, and the red solid lines correspond to the best-fit models, which include both the 
effects of the instrumental variation with time and image position and the planetary transit 
(see §3.2p . In the lower panel for each candidate, we display these data after correcting for 
the instrumental model, normalizing, and binning by 20 minutes. In the case of Kepler-20d 
we co-added the two light curves. The best-fit model for the transits are plotted in red and 
the transit shapes expected from the parameters estimated from the Kepler observations 
(Table [2]) are over-plotted as dashed green lines. The transit depths measured in the Spitzer 
and Kepler bandpasses agree to better than 1 a. 
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3.3. Spectroscopy 

We obtained 30 high resolution spectra of Kepler-20 between UT 2009 Au gust 30 and 



2011 June 16 using the HIRES spectrometer on the Keck I 10-m telescope (jVogt et al. 



1994 ). We took spectra with the same spectrometer set up of HIRES, and with the same 



spectroscopic analysi s, that we normally use for precise Doppler work of nearby FGK stars 



( jjohnson et al.l 1201 ll ). which typically yields a Doppler precision of 1.5 ms~^ for slowly 
rotating FGKM stars. Typical exposure times ranged from 30—45 minutes, yielding an 
SNR of 120 per pixel (1.3 kms~^). The first 9 observations were made with the B5 decker 
(0'.'87 X 3'.'0) that does not permit moonlight subtraction. The remaining 21 observations 
were made with the C2 decker (0'.'87 x 14'.'0) that permits sky subtraction. The internal 
errors were estimated to be between 1.5—2 ms~^. We augmented these uncertainties by 
adding a jitter term of 2.0 ms~^ in quadrature. The earlier 9 observations are vulnerable to 
modest contamination from moonlight, and we have further augmented the uncertainties for 
these 9 values by adding in quadrature a term of 2.7 ms~^, which is based on the ensemble 
performance of stars similarly affected for this magnitude. The final uncertainties range from 
2.5—4 ms~^. The estimated RVs and uncertainties are given in Table [51 We also undertook 
a study of these spectra to determine the spectral line bisectors with the goal of placing 
limits on these sufficient to preclude astrophysical false positives. However, we found that 
the scatter in the bisector centers was somewhat larger than the RV variations, rendering 
the RV detection, while sufficient for mass constraint, inconclusive for confirmation. We 
therefore undertook the statistical study described in §4). 

We obtained two additional spectra of Kepler-20 with HIRES on UT 2009 September 
08 and 2010 August 24. These observations were gathered without the iodine cell, so that 
the extracted stellar spectrum could be used for the template for the RV analysis and for 
estimating stellar parameters. The exposure time for the first spectrum was 30 minutes 
yielding a typical SNR of 90 and the Bl decker (i? = 60, 000) was used. The observing con- 
ditions were slightly better when the second spectrum was obtained and a one hour exposure 
yielded a SNR of 140. The taller B3 decker {R = 60, 000) was used to carry out a better 
sky subtraction. We analyzed these two spectra using the LTE spectral synthesis analysis 



software Spectroscopy Made Easy (SME; IValenti &: Piskunovlll996 



Valenti fc Fischei]l2005f ) 



to estimate the values of Tgff , log g [Fe/H] and v sin We found that the estimates from each 
spectrum were consistent to within 1 a, and hence we averaged our two estimates to obtain 
Teff = 5455 ± 44 K, log^ = 4.4 ± 0.1, [Fe/H]= 0.01 ± 0.04, and wsinz < 2 kms'^ the errors 
listed are those resulting from the analysis of each individual spectrum, and we have refrained 
from assuming a decrease by a factor of a/2. We also proceeded to measure the fiux in the 
cores of the Ca II H and K lines to evaluate the chromospheric activity. We measured that the 
ratio of emission in these lines to the bolometric emission was log-R^K= —4.93 ± 0.05. This 
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estim ate was derived from a Mt. Wilson-style S- value of 0.183 ± 0.005 fllsaacson fc Fischer 
20101), using the measured color B — V = 0.725. The logfZ^j^ value suggests a low activity 
level for a star of this s pectral type, which is co nsistent with the measured y sin i < 2 km s'^. 
Using the relations of iNoyes et al.l (119841 ) and iMamajek &: Hillenbrandl (120081 ) . we infer a 
rotation period of 31 days. 

We also gathered three moderate signal-to-noise ratio, high-resolution spectra of Kepler- 
20 for reconnaissance purposes, two with the Fibre-fed Ech elle Spectrograph (^FIES) a t the 
2.5 m Nordic Optical Telescope (NOT) at La Palma, Spain (jDjupvik fc Andersenll2010l ) and 
one with the Tull Coude Spectrograph on the McDonald observatory 2.7m Harlan Smith 
Telescope. The FIES spectra were taken on 2009 August 5 and 6 using the medium and high 
resolution fibers resulting in a resolution of 46,000 and 67,000, respectively. Each spectrum 
has a wavelength coverage of approximately 360 — 740 nm. The McDonald spectrum was 
taken on 2010 October 25, with a spectral resolution of 60,000. This spectrum was exposed 
to a SNR of 55 per resolution element for the specific purpose of deriving reliable atmospheric 
parameters for the star. 

As an independent check on the parameters derived from the SME an alysis of the 



Keck/ HIRE S data describe d abov e, we derived stellar parameters following [Torres et al. 
(I2OO2I ) and iBuchhave et al.l (120 101 ). As part of this analysis, we employed a new fitting 



scheme that is currently under development by L. Buchhave, allowing us to extract precise 
stellar parameters from the spectra. We analyzed the two FIES spectra, the McDonald 
spectrum and the three HIRES template spectra. These results were found to be consistent 
within the errors. Taking the average of the stellar parameters from the different instruments 
yielded the following parameter estimates: Tes= 5563 ± 50 K, log (7= 4.52 ± 0.10, [m/H] = 
+0.04 ± 0.08, and vsmi= 1.80 ± 0.50 kms~^, in agreement with the parameters from SME 
within the uncertainties. The average systemic radial velocity of the six observations was 
—21.87 ± 0.96 kms~^ on the lAU standard scale, which includes the correction for the 
gravitational redshift of the Sun. 

We note that the two analyses yielded consistent results for log (?, metallicity, and v sini, 
but that the estimates of T^s differed by twice the formal error. Hence we elected to adopt 
the results of the SME analysis for our final values, but we increased the uncertainty on T^s 
to 100 K to reflect the difference between the two estimates. We list our estimates for the 
spectroscopically determined parameters in Table [H 
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Table 5. Keck HIRES Radial Velocity Measurements for Kepler-20 



Date of Obs. [BJD] 


Radial Velocity [ms ^] 


Uncertainty [ms ^J'* 


z4ooU7o.oo57io 


—5.86 


3.78 


2455U74.»Dllciy 


—2.94 


3.74 


2455U75.9UDD7e 


1.09 


3.75 


24o5U7o.8oo79z 


2.57 


3.84 


z4ooU/ 7.yU ( Oo4 


11.08 


3.82 


z4ooUoi.yo2Do7 


— 1.51 


3.83 


24o5Uo2.oci2oo7 


—9.53 


3.78 


2455083.888144 


-7.17 


3.84 


2455084.893633 


4.59 


3.88 


2455134.766132 


-14.15 


3.88 


2455314.027873 


2.05 


2.54 


2455320.085390 


-9.47 


2.59 


2455321.072318 


-4.55 


2.56 


2455345.011208 


-8.71 


2.82 


2455351.073218 


0.41 


2.70 


2455352.043110 


-1.84 


2.50 


2455372.898433 


1.71 


2.62 


2455374.967141 


-16.20 


2.62 


2455378.007402 


-2.59 


2.53 


2455380.950103 


-0.62 


2.55 


2455403.055761 


7.59 


2.96 


2455407.903131 


-10.15 


2.64 


2455411.959424 


-13.79 


2.65 


2455413.894286 


-2.60 


2.59 


2455433.818795 


4.81 


2.46 


2455435.841063 


10.04 


2.58 


2455440.792117 


-10.94 


2.50 


2455500.817440 


-2.94 


3.34 


2455522.725846 


-6.56 


2.63 


2455724.055721 


-2.35 


2.64 



''Includes jitter of 2 ms ^ 



- 31 - 



4. Validation of the Planets Kepler-20b, Kepler-20c, and Kepler-20d 

While the analysis of the radial- velocity ( §3.3p data yielded detections for Kepler-20b and 
Kepler-20c, we found that our analyses of the bisector spans were not sufficient to conffim the 
planetary origin of those variations. Moreover, for Kepler-20d, K00070.04, and K00070.05 
there is no Doppler detection. We therefore rely on a fundamentally different technique 
to establish which, if any, of thes e signals can be persuasively attributed to planets. As 



explained in iLissauer et al.l ( 1201 Ibf ) , when dynamical confirmation of a planet candidate by 
the radial velocity method or by transit timing variations cannot be achieved, we attempt 
instead to validate the candidate by tabulating all viable false positives (blends) that could 
mimic the signal. We then assess the likelihood of these blends, and compare it with an a 
priori estimate of the likelihood that the signal is due to a true planet. We consider the 
signal to be validated when the likelihood of a planet exceeds that of a false positive by a 
sufficiently large ratio, typically at least 300 (i.e. 3 a). 

Our tabulation of the viable scenarios resulting from blends w as accomplished with the 
BLENDER algorithm ( Torres et al. 2004, 2011 : Fressin et al. 2011al lbl) combined with some of 



the follow-up observations described earlier (high-resolution imaging, centroid motion analy- 
sis, spectroscopy, and Spitzer observations). BLENDER attempts to ffi the Kepler photometry 
with a vast array of synthetic light curves generated from blend configurations consisting of 
chance alignments with background or foreground eclipsing binaries (EBs), as well as eclips- 
ing binaries physically associated with Kepler-20 (hierarchical triples). We also considered 
cases in which the second star is eclipsed by a larger planet, rather than by another star. A 
wide range of parameters is explored for the eclipsing pair, as well as for the relative distance 
separating it from the target. Scenarios giving poor ffis to the data (specifically, a value 
that indicates a discrepancy of at least 3 a worse than that corresponding to the transiting 
planet model) are considered to be ruled out. For full details of this technique we refer the 
reader to the above sources. 

The combination of the shorter periods and deeper transits for Kepler-20b and Kepler- 
20c results in higher SNR for those signals compared to the others. Consequently the shape 
of the transit is better defined, and this information makes it easier to reject false positives 
with BLENDER, as we show below. The transit depths of K00070.04 and K00070.05 are only 
82 and 101 parts per million; this renders these signals far more challenging to validate, and 
we find below that we are currently not able to demonstrate unambiguously that these two 
signals are planetary in origin. Kepler-20d is similar in depth to Kepler-20b, but due to its 
longer orbital period, far fewer transits have been observed. This results in a lower SNR in 
the phase light curve. We begin by describing this case. 

Figure M illustrates the BLENDER results for Kepler-20d. The three panels represent 
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cuts through the space of parameters for blends consisting of background EBs, background 
or foreground stars transited by larger planets, and physically associated triples. In the 
latter case we find that the only scenarios able to mimic the signal are those in which 
the companion star is orbited by a larger planet, rather than another star. The orange- 
red-brown-black shaded regions correspond to different levels of the difference between 
blend models and the best transiting planet fit to the Kepler data, expressed in terms of 
the statistical significance of the difference (a). The 3 a level is represented by the white 
contour, and only blends inside it (< 3cr) are considered to give acceptable fits to the 
Kepler photometry. Other constraints further restrict the area allowed for blends. The 
green hatched areas are excluded because the EB is within one magnitude of the target in 
the Kp band, and would generally have been noticed in our spectroscopic observations. The 
blue hatched areas are also excluded because the overall color of the blend is either too red 
(left in the top two panels) or too blue (right) compared to the mea sured Sloan-2MAS S 
r — Ks color of Kepler-20, as listed in the Kepler Input Catalog (KIC; iBrown et al.ll201ll ). 
Additionally, Spitzer observations rule out blends involving EBs (or star-|-planet pairs) with 
stars less massive than about 0.78 Mq (gray shaded area to the left of the vertical dotted line), 
because the predicted depth of the transits in the 4.5 fim bandpass of Warm Spitzer would 
be more than 3 a larger than our Spitzer observations indicate. Note that the combination 
of these constraints rules out all physically associated triple configurations for Kepler-20d, 
so that only certain blend scenarios involving background EBs or background/foreground 
stars transited by larger planets present suitable alternatives to a true planet model. 



We estimate the frequency of these remaining blends following [Torres et al.l ( 120111 ) and 
Fressin et al.l ( l2011al ). as the product of three factors: the expected number density of stars in 
the vicinity of Kepler-20, the area around the target within which we would miss such stars, 
and an estimate of how often we expect those stars to be in EBs or be transited by a larger 
planet of the right characteristics (specified by the stellar masses, planetary sizes, orbital 
eccentricities, and other characteristics as tabulated by B LENDER). For t he nu mber densities 
we appeal to the Besangon Galactic structure model of iRobin et al.l ( 120031 ). Constraints 
from our high-resolution imaging (see §3.1.ip allow us to estimate the maximum angular 
separation (pmax) at which blended stars would be undetected, as a function of brightness. 
We derive our estimates of the frequencies of E Bs and larger transiting planets iriyolve d in 
blends from recent studies by the Kepler Team (jSlawson et al.l 1201 It iBorucki et al.ll201l[ ). in 
the same way as done for our earlier studies of Kepler-9 d, K epler-10 c, and Kepler-llg (see 
Torres et al.l 1201 It iFressin et al.l 1201 lat iLissauer et al.ll2011bl ). 



The results of our calculations for Kepler-20d, performed in half-magnitude bins, are 
shown in Table |6] separately for background EBs and for background or foreground stars 
transited by a larger planet. The first two columns give the Kp magnitude range of each bin. 
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and the magnitude difference ^Kp relative to the target, calculated at the faint end of each 
bin. Column 3 reports the stellar density near the target, subject to the mass constraints from 
BLENDER as shown in Figure O Column 4 gives the maximum angular separation at which 
background stars would escape detection in our imaging observations. In this particular case 
those observations are more constraining than the 3 a exclusion limit set by our analysis 
of the flux centroids (0'.'65; see §3.1.11) . The product of the area imphed by pmax and the 
densities in the previous column are listed in column 5, in units of 10~^. Column 6 is the 
result of multiplying this number of stars by the frequency of suitable EBs (/eb = 0.78%; see 



Fressin et al.ll2011a] ). A similar calculation is performed for blends involving stars transited 
by larger planets, and is presented in columns 7-10, using /pianet = 0.18%. The latter is the 
frequency of planets in th e radius range allow ed by BLENDER for these types of scenarios, 



which is 0.4-2.0 i?jup (see iBorucki et al.ll201ll ). The sum of the contributions in each bin 



is given at the bottom of columns 6 and 10. The total number of blends (i.e., the blend 
frequency) we expect a priori is reported in the last line of the table by adding these two 
numbers together, and is approximately BF = 6.0 x 10"''. 

We now compare this estimate with the likelihood that Kepler-20d is a true transiting 
planet (planet prior). T o calculate the planet prior we again make use of the census of 1,235 



candidates reported by iBorucki et al.l (120111 ) among the 156,453 Kepler targets observed 
during the first four months of operation of the MissioiJ^. We count 100 candidates that are 
within 3 cr of the measured radius ratio of Kepler-20d, implying an a priori transiting planet 
frequency of PF = 100/156,453 = 6.4 x 10~^. The likelihood of a planet is therefore several 
orders of magnitude larger than the likelihood of a false positive (PF/BF = 6.4 x 10~^/6.0 x 
10"'' ~ 1100), and we consider Kepler-20d to be validated as a planet with a high degree of 
confidence. 

The transit signals from Kepler-20b and Kepler-20c are better defined, and as a result 
BLENDER is able to rule out all scenarios involving background EBs consisting of two stars, 
as well as all physically associated triples. This reduces the blend frequencies by one to two 
orders of magnitude compared to Kepler-20d. For Kepler-20c, Spitzer observations are avail- 
able as well, although the constraints they provide are redundant with color information also 
available for the star, which already rules out contaminants of late spectral type. The areas 
of parameter space in which BLENDER finds false positives providing acceptable fits to the 
photometry are shown in Figure [TOl The detailed calculations of the blend frequencies for 



^While these 1,235 candidates have no t yet been confirmed as tr ue planets, the rate of false positives is 
expected to be quite low (10% or less; see lMorton fc Johnsonll2011l ). so our results will not be significantly 
af fected by the assump tion that all of the candidates are planets. We further assume here that the census 
of IBorucki et al.l (j201l[) is complete at these planetary radii. 
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Kepler-20b and Kepler-20c are presented in Table [7] and Table [H respectively, using appro- 
priate ranges for the larger planets orbiting the blended stars as allowed by BLENDER, along 
with the corresponding transiting planet frequencies specified in the headings of column 6. 

Planet pr i ors fo r these two candidates were computed as before using the catalog of 



Borucki et al.l (120 111 ). We count 52 cases in that list within 3 cr of the measured radius ratio 
of Kepler-20b, leading to an a priori planet frequency of 52/156,453 = 3.3 x 10~^. This is 
nearly 20,000 times larger than the blend frequency given in Table [7] (BF = 1.7 x 10"*^). For 
Kepler-20c the planet prior based on the measured radius ratio is 28/156,453 = 1.8 x 10""^, 
which is approximately 10^ times larger than the likelihood of a blend. Therefore, both 
Kepler-20b and Kepler-20c are validated as planets with a very high degree of confidence. 

We carried out similar calculations for the candidates K00070.04 and K00070.05. The 
transit signals of these two candidates are much more shallow than those of Kepler-20b, 
Kepler-20c, and Kepler-20d. As a result, the constraint on the shape of the transit is con- 
siderably weaker than in the cases described above, and many more false positives than 
before are found with BLENDER that provide acceptable fits within 3 a of the quality of a 
planet model. Additionally, neither of these candidates were observed with Spitzer, so the 
constraint on the near-infrared depth of the transit that allowed us to rule out some of the 
blends for Kepler-20d is not available here. In particular, physically associated stars tran- 
sited by a larger planet cannot all be ruled out, and this ends up contributing significantly 
to the overall blend frequency. We conclude that the BLENDER methodology as implemented 
above is insufficient to validate either K00070.04 or K00070.05, and we defer this issue to a 



subsequent study (IFressin et al.ll2012l ). 



Table 6. Blend frequency estimate for Kepler-20d (K00070.03). 



Blends Involving Stellar Tertiaries Blends Involving Planetary Tertiaries 



Kp Range 




Stellar Density^ 


pmax 


Stars 


EBs 


Stellar Density^ 


pmax 


Stars 


Blends (xlO 


(mag) 


(mag) 


(per sq. deg) 


(") 


(xlO~^) 


/eb = 0.78% 


(per sq. deg) 


(") 


(x 10^^) 


Rp e [0.4-2.0 iJjup], 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


12.5-13.0 


0.5 


















13.0-13.5 


1.0 


















13.5-14.0 


1.5 


28 


0.075 


0.038 


0.0003 


273 


0.075 


0.372 


0.0007 


14.0-14.5 


2.0 


77 


0.093 


0.164 


0.0013 


423 


0.093 


0.901 


0.0016 


14.5-15.0 


2.5 


119 


0.11 


0.349 


0.0028 


572 


0.11 


1.678 


0.0024 


15.0-15.5 


3.0 


238 


0.13 


0.975 


0.0077 


897 


0.13 


3.675 


0.0053 


15.5-16.0 


3.5 


532 


0.15 


2.902 


0.0229 


1183 


0.15 


6.452 


0.0093 


16.0-16.5 


4.0 


1321 


0.20 


12.81 


0.1012 


1675 


0.20 


16.24 


0.0234 


16.5-17.0 


4.5 


1593 


0.25 


24.13 


0.1907 










17.0-17.5 


5.0 


1295 


0.30 


28.25 


0.2232 










17.5-18.0 


5.5 


















18.5-19.0 


6.0 


















19.0-19.5 


6.5 


















19.5-20.0 


7.0 


















20.0-20.5 


7.5 


















20.5-21.0 


8.0 


















Totals 




5203 




69.62 


0.5511 


5023 




29.32 


0.0526 



Blend frequency (BF) = (0.5511 + 0.0526) X lO^^ ~ 6.04 X 10"^ 



''The number densities in Columns 3 and 7 differ because of the different secondary mass ranges permitted by BLENDER for the two kinds of blend scenarios, as shown 
in the top two panels of Figures [9] 



Note. — 



Magnitude bins with no entries correspond to brightness ranges in which all blends are ruled out by a combination of BLENDER and other constraints. 
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Table 7. Blend frequency estimate for Kepler-20b (K00070.02). 



Blends Involving Planetary Tertiaries 



Kp Range 


AKp 


Stellar Density 


pmax 


Stars 


Blends (xlO 


(mag) 


(mag) 


(per sq. deg) 


(") 


(xlO-6) 


Rp G [0.27- 1.81 Rjup], 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


12.5-13.0 


0.5 










13.0-13.5 


1.0 










13.5-14.0 


1.5 


271 


0.075 


0.370 


0.00068 


14.0-14.5 


2.0 


379 


0.093 


0.807 


0.00149 


14.5-15.0 


2.5 


404 


0.11 


1.185 


0.00218 


15.0-15.5 


3.0 


464 


0.13 


1.901 


0.00351 


15.5-16.0 


3.5 


498 


0.15 


2.716 


0.00501 


16.0-16.5 


4.0 


255 


0.20 


2.473 


0.00456 


16.5-17.0 


4.5 










17.0-17.5 


5.0 










17.5-18.0 


5.5 










18.0-18.5 


6.0 










18.5-19.0 


6.5 










19.0-19.5 


7.0 










19.5-20.0 


7.5 










20.0-20.5 


8.0 










Totals 




2669 




9.452 


0.0174 



Blend frequency (BF) = 1.74 x 10" 



Note. — Magnitude bins with no entries correspond to brightness ranges in which all blends are ruled out 
by a combination of BLENDER and other constraints. 



Table 8. Blend frequency estimate for Kepler-20c (K00070.01). 



Blends Involving Planetary Tertiaries 



Kp Range l^Kp Stellar Density pmax Stars Blends (xlO^^') 

(mag) (mag) (per sq. deg) (") (xlO-6) ijj, e [0.39-1.95 Rj^p] , /pi3,„ = 0.18% 

(1) (2) (3) (4) (5) (6) 



12.5-13.0 


0.5 










13.0-13.5 


1.0 










13.5-14.0 


1.5 


221 


0.075 


0.300 


0.00056 


14.0-14.5 


2.0 


274 


0.093 


0.551 


0.00108 


14.5-15.0 


2.5 










15.0-15.5 


3.0 










15.5-16.0 


3.5 










16.0-16.5 


4.0 










16.5-17.0 


4.5 










17.0-17.5 


5.0 










17.5-18.0 


5.5 










18.0-18.5 


6.0 










18.5-19.0 


6.5 










19.0-19.5 


7.0 










19.5-20.0 


7.5 










Totals 




495 




0.851 


0.00164 



Blend frequency (BF) = 1.64 x lO"' 



Note. — Magnitude bins with no entries correspond to brightness ranges in which all blends are ruled out 
by a combination of BLENDER and other constraints. 
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0.2 0.4 0.6 0.8 

Secondory moss (Mg) 



Fig. 9. — BLENDER goodness-of-fit contours corresponding to three different blend scenarios 
for Kepler-20d: background EBs {top left), background or foreground stars transited by a 
larger planet {top right), and stars physically associated with the target that are transited 
by a larger planet {bottom). The mass of the intruding star (referred to in the BLENDER 
nomenclature as the secondary) is shown along the horizontal axis, and the distance between 
this star and the target is shown on the vertical axis of the top two panels, expressed for 
convenience in terms of the distance modulus difference. The vertical axis in the bottom 
panel shows the sizes of the planets (tertiaries) orbiting physically associated stars. Viable 
blend models are those giving fits with values within 3 a of the best planet fit, and lie 
inside the white contours. Other colored areas outside the white contours indicate regions 
of parameter space with increasingly worse fits to the data. Blends excluded by our Spitzer 
constraints are shown with the shaded gray area (secondary masses < 0.78 M0). Blue cross- 
hatched areas indicate regions in which blends are excluded because they are either too red 
(left) or too blue (right) compared to the measured r — Kg color of Kepler-20. Blend scenarios 
in the green cross-hatched areas are also ruled out because they are within AKp = 1.0 of 
the brightness of the target, and would have been detected spectroscopically. The dashed 
diagonal green lines in the top panels mark the faintest blends that give acceptable fits to 
the light curve, corresponding to AKp ^ 5 in the top left panel and AKp ^ 4 in the top 
right panel. As a result of the combined constraints from our Spitzer observations, color 
index, and brightness, all physically associated triples are excluded. 
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Fig. 10. — BLENDER constraints for Kepler-20b (top panels) and Kepler-20c (bottom), show- 
ing chance ahgnments with a star+planet pair on the left and physically associated stellar 
companions transited by a larger planet on the right. See Figure [9] for the meaning of the 
various lines. The space of parameters for background EBs is not shown as all of those 
scenarios provide very poor fits to the transit light curve, and are ruled out. All blends 
involving physically associated stars transited by a larger planet ( right panels ) are excluded 
by a combination of spectroscopic constraints (specifically, on the absence of a secondary 
spectrum) and color constraints. 
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5. Constraints on Transit Times and Long-term Stability 



In this section we discuss the transit times and long-term stability of the system of 
planets orbiting Kepler-20. Both are consi stent with the planet in t erpretation for all 5 
candidates: transit timing variations (TTVs; iHolman fc Murrayl 120051 : lAgol et al.ll2005l ) are 
not seen or expected, and the system is expected to be stable over long timescales. 

The individual transit times are measure d by allowing a t emplate transit light curve to 
slide in time to fit the data for each transit (IFord et al.ll201ll ). The resulting transit times 
are given in Table O Aside from slightly more scatter than expected from the formal error 
bars there is no indication of perturbations such as coherent pat terns. Such exces s scatter is 
not atypical of transit times measured by the standard pipeline (jPord et al. 2011 ). Thus, we 
find no evidence for dynamical interactions among either the transiting planets or additional, 
non-transiting planets. 

To calculate predicted transit times, we numerically integrate our baseline model, which 
consists of a central star of mass 0.912 M© surrounded by planets with periods and epochs 
given in Table |2] (given at dynamical epoch BJD 2454170), and with masses of 8.7, 0.65, 
16.1, 1.1 and 8.0 M© (from least to greatest orbital period), corresponding to the best-fit 



masses for Kepler-20b and Kepler-20c, a guess of Mp = M© (_Rp/_R©)^ °^ (ILissauer et al. 



2011bf ) for Kepler-20d, and masses giving Earth's density for K00070.04 and K00070.05. 
(We remind the reader that we have not, in this paper, validated these two candidates as 
planets. However, considering them as such for the purposes of evaluating dynamical stability 
is the conservative choice, since t he presence of 5 pla nets, as opposed to 3, is more likely 
to induce dynamical instabilities. iFressin et al.l (|2012[ ) presents the validation of K00070.04 
and K00070.05, gives their sizes, from which the masses above are derived, and discusses 
constraints on their masses.) The orbits are chosen to be initially circular, coplanar, and 
edge-on to the line of sight. The root mean square deviations of the model transit times from 
the best-fit linear ephemeris projected over 8 years are approximately 3s, 76s, 9s, 95s and 
10s (from least to greatest orbital period), all significantly smaller than the measurement 
precision shown in Figure [TTl 

Next, we investigated long-term stabilit y for this system by integrating the baseline 
model with the hybrid algorithm in Mercury ( lChamberslll999l ). As no close encounters wer e 
recorded, this algorithm reduced to the symplectic algorithm (jWisdom fc HolmanI Il99ll ). 
with time steps of 0.1 days, roughly 2.7% of the period of the innermost planet. Over the 
10 Myr integration duration, there were no indications of instability. The orbital eccentrici- 
ties fiuctuated on the scale between approximately 3 x 10~^ (for Kepler-20d) and 0.001 (for 
K00070.04). We conclude that plausible, low-eccentricity models for the system are stable 
over long timescales. 
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Finally, we performed an ensemble of N-bo dy integrations using the time-symmetric 4th 
order Hermite integrator ( Kokubo et al. 1998 ) implemented in Swarm-N(if| to estimate the 
maximum plausible eccentricity for each planet consistent with long-term stability. For each 
N-body integration, we set four planets on circular orbits and assigned one planet a non- 
zero eccentricity. The eccentricity and pericenter directions for the planet on a non-circular 
orbit were drawn from uniform distributions. The maximum for the uniform distribution 
of eccentricities was chosen to be slightly larger than necessary for its orbit to cross one of 
its neighbors. We report Cmax; the maximum initial eccentricity that resulted in a system 
with no close encounters (within one mutual Hill radius) and semi-major axes (in a Jacobi 
frame) which varied by less than 1% for the duration of the integrations. Based on 100 
integrations per planet and relatively short integrations (10^ years), we estimate Cmax to be 
0.19, 0.16, 0.16, 0.38 and 0.55 (from smallest to largest orbital period). Technically, we can 
not completely exclude larger eccentricities, due to various assumptions (such as the planet 
masses, coplanarity, prograde orbits, absence of false positives, and the potential for small 
islands of stability at higher eccentricity). Nevertheless, the N-body integrations support 
the assumption of non-crossing orbits, as the vast majority of systems with an eccentricity 
larger than Cmax are dynamical unstable. 



www. astro. ufl.edu/~eford/code/swarm/ 
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Table 9. Transit Times for Kepler-20 We request that this table be published in electronic 

format. 
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.0028 


Kepler-20b 


31 


182.0814 


0.0013 





.0038 


Kepler-20b 


32 


185.7734 


-0.0027 


0, 


.0038 


Kepler-20b 


33 


189.4547 


-0.0176 


0, 


.0103 


Kcplcr-20b 


35 


196.8662 


0.0016 





.0038 


Kcplcr-20b 


36 


200.5665 


0.0059 





.0036 


Kcplcr-20b 


37 


204.2595 


0.0027 





.0026 


Kepler-20b 


38 


207.9564 


0.0035 


0, 


.0032 


Kepler-20b 


39 


211.6447 


-0.0043 





.0030 


Kepler-20b 


40 


215.3402 


-0.0049 


0, 


.0054 


Kepler-20b 


41 


219.0422 


0.0009 


0, 


.0040 


Kcplcr-20b 


42 


222.7391 


0.0017 





.0048 


Kcplcr-20b 


43 


226.4399 


0.0064 





.0052 


Kcplcr-20b 


44 


230.1275 


-0.0022 





.0024 


Kepler-20b 


45 


233.8232 


-0.0025 


0, 


.0031 


Kepler-20b 


46 


237.5259 


0.0041 


0, 


.0043 


Kepler-20b 


47 


241.2178 


-0.0002 


0, 


.0037 
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n 


tn 

BJD-2454900 


-I J- Vn 
W 




w 


Kepler-20b 


48 


OA A 


noi o 


0.0071 





.0035 


Kcplcr-20b 


49 


nAQ 


.DU ( □ 


-0.0028 





.0033 


Kcplcr-20b 


50 




.oUo / 


-0.0007 





.0045 


Kepler-20b 


52 




.Doyo 


-0.0093 


0, 


.0043 


Kepler-20b 


53 


ZOO. 


.4U1 ( 


0.0070 





.0046 


Kepler-20b 


54 


Zui . 


no/1 A 
.uy44 


0.0036 





.0039 


Kepler-20b 


55 


z (v. 


. (iSoZ 


-0.0038 





.0051 


Kcplcr-20b 


56 


Z i^. 


AQ^ A 


-0.0017 





.0042 


Kcplcr-20b 


57 


Z to. 


. ioiO 


0.0024 





.0049 


Kcplcr-20b 


58 


OCT 

Zol. 


.Of Zo 


-0.0026 





.0041 


Kepler-20b 


59 




.O ( o4 


0.0020 


0, 


.0044 


Kepler-20b 


60 


Zoi). 


.ZDoO 


-0.0041 





.0041 


Kepler-20b 


61 


ZtyZ. 


QdKA 
.yD04 


0.0017 





.0053 


Kepler-20b 


62 




.Dooy 


-0.0040 





.0028 


Kcplcr-20b 


63 


oUU. 


.0044 


-0.0015 





.0039 


Kcplcr-20b 


64 




.U04U 


0.0019 





.0035 


Kcplcr-20b 


65 


OU / . 


7A7Q 


-0.0002 





.0040 


Kepler-20b 


66 


oil. 


AAAQ 
.444y 


0.0005 





.0045 


Kepler-20b 


67 


olO. 


.14 / 4 


0.0070 





.0042 


Kepler-20b 


68 


oio. 


.OO40 


-0.0020 


0, 


.0045 


Kepler-20b 


69 


ozz, 


.Do i 4 


0.0048 





.0040 


Kepler-20b 


70 


OZO. 


.ZZ 1 4 


-0.0014 





.0039 


Kepler-20b 


74 




.uiyo 


0.0062 


0, 


.0040 


Kepler-20b 


75 


044, 


. f IDO 


0.0072 
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Kepler-20b 


76 
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.0048 


Kepler-20b 


77 
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Keplcr-20b 


78 


ooo. 
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0.0047 
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Kepler-20b 


79 


359, 


.4904 


-0.0035 





.0043 


Kepler-20b 


80 


363, 


.1916 


0.0016 





.0045 


Kepler-20b 


81 


366, 


.8831 


-0.0031 


0, 


.0041 


Kepler-20b 


82 


370, 


.5829 


0.0006 





.0037 


Kepler-20b 


83 


374, 


.2740 


-0.0044 





.0027 


Kepler-20b 


84 


377, 


.9727 


-0.0018 





.0040 


Kcplcr-20b 


85 


381, 


.6793 


0.0086 





.0039 


Kcplcr-20b 


86 


385, 


.3678 


0.0011 





.0033 


Kcplcr-20b 


87 


389, 


.0654 


0.0025 





.0041 


Kepler-20b 


88 


392, 


.7595 


0.0005 


0, 


.0045 


Kepler-20b 


89 


396, 


.4618 


0.0067 





.0042 


Kepler-20b 


90 


400, 


.1460 


-0.0053 





.0042 


Kepler-20b 


91 


403, 


.8471 


-0.0003 





.0043 


Kcplcr-20b 


92 


407, 


.5457 


0.0022 





.0027 


Kepler-20b 


93 


411, 


.2340 


-0.0056 


0, 


.0049 


Kepler-20b 


94 


414, 


.9291 


-0.0066 





.0036 


Kepler-20b 


95 


418, 


.6325 


0.0007 





.0034 


Kepler-20b 


96 


422, 


.3282 


0.0003 





.0029 
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Kepler- 
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Kepler- 
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.0035 


Kepler- 
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0.0129 
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.0053 


Kepler- 
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.oy 1 / 


-0.0013 
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Kepler- 
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4oo. 
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A QO 


.00o4 


-0.0009 


0, 


.0035 


Kepler- 


■20b 


116 


ACkfi 
490. 




0.0022 


0, 


.0033 


Kepler- 


■20b 


118 


OUo. 


fiA 1 1 
.0411 


-0.0015 


0, 


.0039 


Kepler- 


-20b 


119 


f^n7 
OU ( . 


,OZlO 


-0.0173 


0, 


.0045 


Kcplcr- 


-20b 


120 


Oil. 


.Uo i 1 


0.0022 
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.0055 


Kepler- 


-20b 


121 


K^ A 

014. 


. /olo 


0.0008 


0, 


.0049 


Kepler- 


-20b 
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Oio. 


.4zOO 


-0.0016 


0, 
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Kepler- 
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1 1 Qc^ 

.iiyo 


-0.0037 


0, 


.0043 


Kepler- 


■20b 
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ozo. 


m fin 


-0.0034 


0, 


.0030 


Keplcr- 


-20b 


125 


ozy , 




-0.0001 


0, 


.0044 


Kepler- 


-20b 


126 


533, 


.2147 


0.0031 


0, 


.0046 


Kepler- 


-20b 


127 


536, 


.9038 


-0.0039 


0, 


.0029 


Kepler- 


■20b 


128 


540, 


.6034 


-0.0004 


0, 


.0037 


Kepler- 


■20b 


129 


544, 


.3084 


0.0084 


0, 


.0042 


Kepler- 


■20b 


130 


548, 


.0018 


0.0057 


0, 


.0041 


Kepler- 


■20b 


131 


551, 


.6946 


0.0024 


0, 


.0035 


Kepler- 


-20b 


132 


555, 


.3912 


0.0028 


0, 


,0030 


Kepler- 


-20b 


133 


559, 


.0861 


0.0016 


0, 


.0037 


Kepler- 


-20b 


135 


566, 


.4817 


0.0050 


0, 


.0043 


Kepler- 


■20b 


136 


570, 


.1738 


0.0009 


0, 


.0054 


Kepler- 


■20b 


137 


573, 


.8736 


0.0046 


0, 


.0032 


Kepler- 


■20b 


138 


577, 


.5612 


-0.0039 


0, 


.0040 


Kepler- 


■20b 


139 


581, 


.2539 


-0.0074 


0, 


.0051 


Kepler- 


■20b 


140 


584, 


.9482 


-0.0092 


0, 


.0040 


Kepler- 


-20b 


141 


588, 


.6502 


-0.0032 


0, 


.0035 


Kepler- 


-20b 


142 


592, 


.3510 


0.0014 


0, 


.0050 


Kepler- 


■20b 


143 


596, 


.0490 


0.0033 


0, 


.0034 


Kepler- 


■20b 


144 


599, 


.7493 


0.0075 


0, 


.0044 
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ID 


n 


in 


TTV„ 








BJD-2454900 


(d) 


(d) 



Kepler-20b 


145 


603.4392 


0.0013 


0, 


,0056 


Kcplcr-20b 


146 


607.1344 


0.0004 


0, 


.0038 


Kcplcr-20b 


147 


610.8329 


0.0027 


0, 


.0043 


Kepler-20b 


148 


614.5286 


0.0023 


0, 


.0032 


Kepler-20b 


149 


618.2205 


-0.0019 


0, 


.0027 


Kepler-20b 


150 


621.9156 


-0.0029 


0, 


.0031 


Kepler-20b 


151 


625.6180 


0.0033 


0, 


.0039 


Kepler-20b 


152 


629.3115 


0.0007 


0, 


,0039 


Kepler-20b 


153 


633.0118 


0.0049 


0, 


,0050 


Kepler-20b 


154 


636.6958 


-0.0072 


0, 


,0043 


Kepler-20b 


155 


640.3937 


-0.0054 


0, 


.0045 


Kepler-20b 


157 


647.7922 


0.0008 


0, 


.0046 


Kepler-20b 


158 


651.4935 


0.0060 


0, 


.0042 


Kepler-20b 


163 


669.9606 


-0.0075 


0, 


,0041 


Kepler-20b 


164 


673.6626 


-0.0017 


0, 


,0094 


Kepler-20b 


165 


677.3552 


-0.0052 


0, 


,0037 


Kepler-20b 


166 


681.0550 


-0.0015 


0, 


,0062 


Kepler-20b 


167 


684.7541 


0.0014 


0, 


,0030 


Kepler-20b 


168 


688.4526 


0.0038 


0, 


.0046 


Kepler-20b 


169 


692.1463 


0.0014 


0, 


,0017 


Kepler-20b 


171 


699.5324 


-0.0047 


0, 


,0039 


Kcplcr-20b 


172 


703.2330 


-0.0002 


0, 


,0067 


Kepler-20b 


173 


706.9307 


0.0014 


0, 


,0038 


Kepler-20b 


174 


710.6210 


-0.0045 


0, 


.0032 


Kepler-20b 


175 


714.3175 


-0.0041 


0, 


.0037 


Kepler-20b 


176 


718.0073 


-0.0105 


0, 


,0063 


Keplcr-20b 


177 


721.7094 


-0.0044 


0, 


,0042 


Kepler-20b 


178 


725.4097 


-0.0003 


0, 


,0039 


Kepler-20b 


179 


729.1068 


0.0007 


0, 


,0039 


Kepler-20b 


180 


732.8057 


0.0035 


0, 


,0032 


70.04 




68.9336 + n X 6.098493 






70.04 





68.9295 


-0.0041 


0, 


.0124 


70.04 


1 


75.0316 


-0.0005 


0, 


,0149 


70.04 


2 


81.1487 


0.0181 


0, 


,0116 


70.04 


3 


87.1926 


-0.0364 


0, 


,0176 


70.04 


4 


93.2997 


-0.0278 


0, 


,0117 


70.04 


6 


105.5245 


-0.0001 


0, 


,0103 


70.04 


7 


111.6445 


0.0215 


0, 


,0150 


70.04 


8 


117.6964 


-0.0251 


0, 


.0150 


70.04 


9 


123.8323 


0.0123 


0, 


,0138 


70.04 


10 


129.9095 


-0.0090 


0, 


,0122 


70.04 


11 


135.9953 


-0.0217 


0, 


,0525 


70.04 


12 


142.0980 


-0.0175 


0, 


.0176 


70.04 


13 


148.2067 


-0.0073 


0, 


.0211 


70.04 


14 


154.3205 


0.0080 


0, 


,0156 
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ID 


n 




TTV„ 








BJD-2454900 


(d) 


(d) 



70, 


.04 


15 


160, 


.4350 


0.0240 


0, 


.0122 


70, 


.04 


16 


166, 


.5186 


0.0091 


0, 


.0152 


70, 


.04 


17 


172, 


.6143 


0.0063 


0, 


.0123 


70, 


.04 


18 


178, 


.7111 


0.0047 


0, 


.0141 


70, 


.04 


19 


184. 


.8108 


0.0059 


0, 


.0503 


70, 


.04 


20 


190, 


.9342 


0.0307 


0, 


.0171 


70 


.04 


21 


196, 


.9966 


-0.0054 


0, 


.0137 


70, 


.04 


22 


203, 


.1011 


0.0007 


0, 


.0148 


70, 


.04 


23 


209, 


.2095 


0.0106 


0, 


.0177 


70, 


.04 


24 


215, 


.3090 


0.0115 


0, 


.0181 


70, 


.04 


25 


221, 


.3969 


0.0010 


0, 


.0163 


70, 


.04 


26 


227. 


.5457 


0.0513 


0, 


.0112 


70, 


.04 


27 


233. 


.6113 


0.0183 


0, 


.0132 


70, 


.04 


28 


239, 


.7084 


0.0170 


0, 


.0226 


70, 


.04 


29 


245, 


.7739 


-0.0160 


0, 


,0192 


70, 


.04 


30 


251, 


.8802 


-0.0082 


0, 


,0147 


70, 


.04 


31 


257, 


.9671 


-0.0198 


0, 


,0147 


70, 


.04 


32 


264, 


.0366 


-0.0488 


0, 


.0175 


70, 


.04 


33 


270, 


.1797 


-0.0042 


0, 


.0314 


70, 


.04 


34 


276, 


.2872 


0.0048 


0, 


.0161 


70, 


.04 


35 


282, 


.3650 


-0.0159 


0, 


.0197 


70, 


.04 


36 


288, 


.4550 


-0.0244 


0, 


,0211 


70, 


.04 


37 


294, 


.5290 


-0.0489 


0, 


.0186 


70, 


.04 


38 


300. 


.6801 


0.0037 


0, 


.0163 


70, 


.04 


39 


306, 


.7537 


-0.0211 


0, 


.0169 


70 


.04 


40 


312, 


.9025 


0.0292 


0, 


.0190 


70, 


.04 


41 


318, 


.9743 


0.0024 


0, 


.0131 


70, 


.04 


42 


325, 


.0517 


-0.0186 


0, 


,0129 


70, 


.04 


45 


343, 


.4285 


0.0627 


0, 


.0183 


70, 


.04 


46 


349, 


.4757 


0.0114 


0, 


.0189 


70, 


.04 


47 


355. 


.5391 


-0.0236 


0, 


.0120 


70, 


.04 


48 


361, 


.6421 


-0.0192 


0. 


.0113 


70, 


.04 


49 


367, 


.7477 


-0.0120 


0, 


.0104 


70, 


.04 


50 


373, 


.8366 


-0.0217 


0, 


,0106 


70, 


.04 


51 


379, 


.9503 


-0.0065 


0, 


,0109 


70, 


.04 


52 


386, 


.0611 


0.0059 


0, 


.0128 


70, 


.04 


53 


392, 


.1451 


-0.0086 


0, 


.0139 


70, 


.04 


54 


398. 


.2748 


0.0226 


0, 


.0204 


70, 


.04 


55 


404. 


.3426 


-0.0082 


0. 


.0139 


70 


.04 


57 


416, 


.5261 


-0.0216 


0, 


.0164 


70, 


.04 


58 


422, 


.6521 


0.0059 


0, 


.0142 


70, 


.04 


59 


428, 


.7752 


0.0305 


0, 


.0155 


70 


.04 


60 


434. 


.8565 


0.0133 


0, 


.0164 


70, 


.04 


61 


440, 


.9316 


-0.0101 


0, 


.0151 


70, 


.04 


62 


447. 


.0528 


0.0127 


0, 


.0172 
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w 
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70.04 


63 


453.1351 


-0.0035 


0.0135 


70.04 


64 


459.2399 


0.0028 


0.0119 


70.04 


65 


465.3240 


-0.0116 


0.0147 


70.04 


67 


477.5344 


0.0017 


0.0120 


70.04 


68 


483.6637 


0.0326 


0.0097 


70.04 


69 


489.7101 


-0.0195 


0.0152 


70.04 


70 


495.8219 


-0.0062 


0.0158 


70.04 


71 


501.9295 


0.0029 


0.0135 


70.04 


72 


507.9607 


-0.0644 


0.0155 


70.04 


73 


514.1161 


-0.0075 


0.0148 


70.04 


74 


520.1992 


-0.0229 


0.0127 


70.04 


75 


526.3201 


-0.0005 


0.0121 


70.04 


76 


532.4389 


0.0198 


0.0133 


70.04 


77 


538.5127 


-0.0049 


0.0143 


70.04 


78 


544.6282 


0.0121 


0.0113 


70.04 


79 


550.7134 


-0.0011 


0.0182 


70.04 


80 


556.8084 


-0.0046 


0.0092 


70.04 


82 


568.9873 


-0.0227 


0.0189 


70.04 


83 


575.1346 


0.0261 


0.0154 


70.04 


84 


581.2328 


0.0258 


0.0157 


70.04 


85 


587.2972 


-0.0083 


0.0147 


70.04 


87 


599.4992 


-0.0033 


0.0139 


70.04 


88 


605.6078 


0.0068 


0.0158 


70.04 


89 


611.6560 


-0.0435 


0.0111 


70.04 


90 


617.8120 


0.0141 


0.0169 


70.04 


92 


630.0025 


0.0075 


0.0151 


70.04 


93 


636.1181 


0.0247 


0.0147 


70.04 


94 


642.1929 


0.0010 


0.0106 


70.04 


95 


648.3022 


0.0118 


0.0171 


70.04 


99 


672.6873 


0.0029 


0.0192 


70.04 


100 


678.7764 


-0.0065 


0.0153 


70.04 


101 


684.8819 


0.0006 


0.0144 


70.04 


102 


690.9771 


-0.0027 


0.0165 


70.04 


103 


697.0837 


0.0053 


0.0166 


70.04 


104 


703.1848 


0.0079 


0.0228 


70.04 


105 


709.3466 


0.0713 


0.0171 


70.04 


106 


715.3852 


0.0113 


0.0108 


70.04 


107 


721.4592 


-0.0132 


0.0155 


70.04 


108 


727.5689 


-0.0019 


0.0169 


70.04 


109 


733.6590 


-0.0104 


0.0103 


Keplcr-20c 




71.6076 + n 


X 10.85409 




Kcplcr-20c 





71.6063 


-0.0013 


0.0018 


Kepler-20c 


1 


82.4610 


-0.0007 


0.0027 


Kepler-20c 


2 


93.3240 


0.0083 


0.0033 


Kepler-20c 


3 


104.1716 


0.0017 


0.0024 
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0, 
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.1400 


-0.0002 





.0018 


Kepler-20c 


11 
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12 
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13 


ZiZ, 


71 07 
. i Iz / 
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14 


zzo, 
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Kepler-20c 


15 




.4zUo 
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0, 
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Kepler-20c 


16 


OA ^ 


07/1/1 

.Z / 441: 
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Kepler-20c 


18 


zoo. 
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0.0086 
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Kepler-20c 


19 


077 

Z I 1 . 


.OOOO 


0.0034 





.0018 


Kcplcr-20c 


20 


Zoo, 


.by 11 


0.0017 





.0019 


Kcplcr-20c 


21 




/I 

,0440 


0.0009 





.0018 


Kepler-20c 


22 


OiU, 


, oyoo 


-0.0013 





.0016 


Kepler-20c 


23 


OZl, 


oc;m 


-0.0016 


0, 


.0018 


Kepler-20c 


25 




.youo 


0.0009 





.0018 


Kepler-20c 


26 


oOo, 


CI oo 

.oizy 


-0.0011 





.0014 


Kepler-20c 


27 


oD4, 


.DO ( O 


-0.0003 


0, 
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Kcplcr-20c 


29 


otto. 


Q7f^ T 
. o / 1 


-0.0011 
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Kepler-20c 


30 


'iCi'7 


o'im 

.ZtjUl 
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0, 
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Kepler-20c 


32 


/tic 


QQQO 

.yoyz 
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Kepler-20c 
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7QOQ 
. / oZo 


-0.0103 
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Kepler-20c 


34 
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.ouuo 
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.0012 
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36 


462, 


.3549 


0.0000 
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Kepler-20c 


37 


473, 


.2071 


-0.0019 





.0015 


Kepler-20c 


38 


484, 


.0631 


0.0000 


0, 


.0017 


Kepler-20c 


39 


494. 


.9177 


0.0006 





.0029 


Kepler-20c 


40 


505, 


.7737 


0.0024 





.0011 


Kepler-20c 


41 


516. 


.6244 


-0.0010 





.0016 


Kcplcr-20c 


42 


527, 


.4780 


-0.0015 





.0013 


Kcplcr-20c 


43 


538, 


.3328 


-0.0008 





.0014 


Kcplcr-20c 


44 


549, 


.1875 


-0.0001 





.0021 


Kepler-20c 


45 


560, 


.0406 


-0.0011 


0, 


.0015 


Kepler-20c 


46 


570. 


.8982 


0.0024 





.0014 


Kepler-20c 


47 


581. 


.7495 


-0.0004 





.0021 


Kepler-20c 


48 


592. 


.6034 


-0.0006 





.0034 


Kepler-20c 


49 


603, 


.4558 


-0.0023 





.0033 


Kepler-20c 


50 


614, 


.3101 


-0.0021 


0, 


.0023 


Kepler-20c 


51 


625. 


.1658 


-0.0005 





.0014 


Kepler-20c 


52 


636. 


.0218 


0.0014 





.0014 


Kepler-20c 


53 


646. 


.8745 


0.0000 





.0016 
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(d) 


(d) 


Kcplcr-20c 


55 


668.5835 


0.0009 


0.0016 


Kcplcr-20c 


56 


679.4363 


-0.0005 


0.0025 


Kcplcr-20c 


57 


690.2898 


-0.0010 


0.0023 


Kepler-20c 


58 


701.1456 


0.0017 


0.0023 


Kepler-20c 


59 


711.9975 


-0.0015 


0.0018 


Kepler-20c 


60 


722.8602 


0.0071 


0.0024 


Kepler-20c 


61 


733.7038 


-0.0034 


0.0014 


70.05 




68.219 + n 


X 19.57706 




70.05 





68.2114 


-0.0076 


0.0111 


70.05 


1 


87.7637 


-0.0323 


0.0155 


70.05 


2 


107.3404 


-0.0327 


0.0149 


70.05 


3 


126.9232 


-0.0270 


0.0161 


70.05 


4 


146.5711 


0.0439 


0.0177 


70.05 


5 


166.0304 


-0.0739 


0.0121 


70.05 


6 


185.6525 


-0.0289 


0.0138 


70.05 


7 


205.2436 


-0.0149 


0.0191 


70.05 


8 


224.8237 


-0.0117 


0.0134 


70.05 


9 


244.4134 


0.0009 


0.0259 


70.05 


10 


264.0047 


0.0151 


0.0129 


70.05 


12 


303.1682 


0.0245 


0.0134 


70.05 


13 


322.7048 


-0.0160 


0.0135 


70.05 


14 


342.2962 


-0.0017 


0.0119 


70.05 


15 


361.8658 


-0.0091 


0.0217 


70.05 


16 


381.4339 


-0.0181 


0.0154 


70.05 


17 


401.0316 


0.0025 


0.0157 


70.05 


18 


A nr\ ^ f\ Aa 

420.5948 


-0.0113 


0.0129 


70.05 


19 


A A r\ ^ a A A 

440.1844 


0.0013 


0.0128 


70.05 


20 


459.7616 


0.0014 


0.0193 


70.05 


21 


479.3348 


-0.0024 


0.0110 


70.05 


22 


498.9270 


0.0127 


0.0138 


70.05 


23 


518.4884 


-0.0030 


0.0113 


70.05 


24 


538.0258 


-0.0427 


0.0118 


70.05 


25 


557.6534 


0.0079 


0.0114 


70.05 


26 


577.2165 


-0.0061 


0.0207 


70.05 


27 


596.8058 


0.0062 


0.0157 


70.05 


28 


616.4009 


0.0242 


0.0148 


70.05 


29 


635.9785 


0.0247 


0.0145 


70.05 


31 


675.1005 


-0.0074 


0.0192 


70.05 


33 


714.2609 


-0.0010 


0.0139 


70.05 


34 


733.8455 


0.0065 


0.0101 


Kcplcr-20d 




97.7271 + n 


X 77.61184 




Kepler-20d 





97.7293 


0.0022 


0.0029 


Kepler-20d 


1 


175.3387 


-0.0002 


0.0022 


Kepler-20d 


2 


252.9482 


-0.0026 


0.0040 


Kepler-20d 


5 


485.7834 


-0.0029 


0.0022 
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641.0100 


0.0001 
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0.0013 
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Fig. 11. — Measured and predicted transit timing for the planets of Kepler-20. Left panels: 
Observed times minus calculated times according to a constant-period model {O — C) are 
plotted as points with error bars, versus transit time. The timing simulations using circular, 
coplanar planets with nominal masses are plotted as lines. Right panels: The simulations 
are shown in more detail (30 x zoom of each panel) to show the timescale and structure of 
variations. 
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6. Constraints on the Planetary Compositions and Formation History 



The Kepler-20 system, harboring muhiple sub-Neptune planets with constrained radii 
and masses, informs our understanding of both models of planet formation and the interior 
structure of planets that straddle the boundary between sub-Neptunes and super-Earths. 
The transit radii measured by Kepler and the planetary masses measured (or bounded) 
by radial velocity observations together constrain the interior compositions of Kepler-20b, 
Kepler-20c, and Kepler-20d, as illustr ated by the mass-radius diagram (Figure [121) . We em- 
ploy planet interior structure models ( Rogers fc Seager 2010 : Rogers et al. 201 ll ) to explore 
the range of plausible planet compositions. The interpretation is challenging because we do 
not yet know if these sub- Neptune planets had a stunt ed formation, or if they formed as 
gas giants and then lost significant mass to evaporation ( iBaraffe et al.ll2004l ). This is partly 
owing to the uncertainties involved in atmospheric escape modeling. 

Notably, both Kepler-20c and Kepler-20d require significant volatile contents to account 
for their low mean densities, and cannot be composed of rocky and iron material alone. The 
volatile material in these planets could take the form of ices (H2O, CH4, NH3) and/or H/He 
gas accreted during planet formation. Outgassing of rocky planets releases an insufficient 
quantity of volatiles (no more than 23% H2O and 3.6% H2 relative to the planet mass) to 
account for Kepler-20 c and could account for Kepl e r-20d only in fine-tuned near-optimal 



outga ssing scenarios (lElkins-Tanton fc Seagerl l2008l : ISchaefer fc Fegleyl l2008l : [Rogers et al. 



2OIII ). For Kepler-20c, ices (likely dominated by H2O) would need to constitute the major- 
ity of its mass, in the absence of a voluminous, though low-mass, envelope of light gases. 
Alternatively, a composition with approximately 1% by mass H/He surrounding an Earth- 
composition refractory interior also matches the observed properties of the planet within 
1 a. Intermediate scenarios, wherein both H/He and higher mean molecular weight volatile 
species from ices contribute to the planet mass, are also possible. For Kepler-20d, the 2 a 
upper limit on the planet density demands at least a few percent H2O by mass, or a few 
tenths of a percent H/He by mass. 

The nature of Kepler-20b's composition is ambiguous: Kepler-20b could be terrestrial 
(with the transit radius defined by a rocky surface), or it could support a significant gas 
envelope (like Kepler-20c and Kepler-20d). In the mass-radius diagram (Figure [T2|) . the 
measured properties of Kepler-20b straddle the pure-silicate composition curve that defines 
a strict upper bound to rocky planet radii. If Kepler-20b is in fact a terrestrial planet 
consisting of an iron core surrounded by a silicate mantle, the 1 a limits on the planet mass 
and radius constrain the iron core to be less than 62% of the planet mass. In particular, 
an Earth-like composition (30% iron core, 70% silicate mantle) is possible and matches the 
observational constraints to within 1 a, but a Mercury-like composition (70% iron core. 
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30% silicate mantle) is not acceptable. Alternatively, Kepler-20b may harbor a substantial 
gas layer like its sibling planets Kepler-20c and Kepler-20d at larger orbital semi-major 
axes, and/or contain a significant component of astrophysical ices such as H2O. The 1 a 
lower limits on the planet density constrain the fraction of Kepler-20b's mass that can be 
contributed by H2O (< 55%) and H/He (< 1%). 

Given their high levels of stellar irradiation (the semi-major axes of all five Kepler-20 
planets are smaller than that of Mercury), atmospheric escape likely played an important 
role sculpting the compositions of the Kepler-20 planets. Planet compositions with low 
mean molecular weight gas envelopes would be especially susc eptible to mass loss. Using a 



model for energy limited escape from hydrogen-rich envelopes (iLecavelier Pes Etangsll2007l ). 
we estimate that Kepler-20b would be losing on the order of 4 x 10^ kgs~^, which corre- 
sponds to 0.02 M^Gyr^^. Following the same approach, the estimated hydrogen mass loss 
rates for Kepler-20c and Kepler-20d are 2 x 10^ kgs'^ (0.01 MeGyr^^) and 8 x 10^ kgs"^ 
(0.0004 M^Gyr^^), respectively. Our theoretical understanding of atmospheric escape from 
highly irradiated super-Earth and sub-Neptune exoplanets is very uncertain, and higher mass 
loss rates are plausible (especially at earlier times when the host star was more active). It is 
intriguing that Kepler-20b, with its shorter orbital period and greater vulnerability to mass 
loss, also has a higher mean density than Kepler-20c and Kepler-20d. More detailed mod- 
eling may constrain Kepler-20b's compositional history and the extent to which its relative 
paucity of volatiles can be attributed to evaporation. 

The Kepler-20 pl anetary system shares several remarkable attributes with Kepler-11 



(ILissauer et al.l l2011bl ). namely the presence of multiple transiting low-density low-mass 
planets in a closely spaced orbital architecture. The Kepler-20 system is less extreme than 
Kepler-11 in the realms of both low planet densities (Figure [T2|) and dynamical compactness 
(the Kepler-11 planets exhibit TTVs while the Kepler-20 planets do not). 

A striking feature of the Kepler-20 planetary system is the presence of Earth-size rocky 
planet candidates interspersed between volatile-rich sub-Neptunes at smaller and larger or- 



bital s emi-major axes, as also seen in Kepler candidate multi-planet systems (ILissauer et al. 



2011al ). Assuming that both K00070.04 and K00070.05 are planets, the distribution of the 
Kepler-20 planets in orbital order is as follows: Kepler-20b (3.7 days, 1.9 -R©), K00070.04 
(6.1 days, 0.9 i?®), Kepler-20c (10.9 days, 3.1 K00070.05 (19.6 days, 1.0 i?e), and 

Kepler-20d (77.6 days, 2.8 -R©). Given the radii and irradiation fluxes of the two Earth-size 
planet candidates, they would not retain gas envelopes. The first, second, and fourth planets 
have high densities indicative of solid planets, while the other two planets have low densities 
requiring significant volatile content. The volatile-rich third planet, Kepler-20c dominates 
the inner part of the Kepler-20 system, by holding much more mass than the other three inner 



planets put together. In the Solar System, the terrestrial planets, gas-giants, and ice giants 
are neatly segregated in regions with increasing distance from the sun. Planet formation 
theories were developed to retrodict these S olar System composition trends (e.g.. 



Safronov 



19691 : IChambersI l20ld : iD'Angelo et al.ll2010l ). In the Kepler-20 system, the locations of the 
low-density sub-Neptunes that are rich in water and/or gas, and the Earth-size planet can- 
didates does not exhibit a clean ordering with orbital period, challenging the conventional 
planet formation paradigm. In situ assembly may form multi-planet systems with close-in 
hot-Neptunes and super-Earths, provided the initial proto planetary disk contained massive 
amounts of solids (~ 50-100 M®) within lAU of the star (IHansen &: Murrayll201ll ). 



Kepler was competitively selected as the tenth Discovery mission. Funding for this 
mission is provided by NASA's Science Mission Directorate. The authors thank many people 
who gave so generously of their time to make this mission a success. This work is also based 
in part on observations made with the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. 
Support for this work was provided by NASA through an award issued by JPL/Caltech. 
We would like to thank the Spitzer staff at IPAC and in particular Nancy Silbermann for 
scheduling the Spitzer observations of this program. Some of the data presented herein 
were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership 
among the California Institute of Technology, the University of California and the National 
Aeronautics and Space Administration. The Observatory was made possible by the generous 
financial support of the W. M. Keck Foundation, (c) 2011 all rights reserved. 
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Fig. 12. — Mass-radius relationships of small transiting planets. The three confirmed planets 
in the Kepler-20 system are highlighted in green. Kepler-20b and Kepler-20c are plotted with 
error bars delimiting the 1 a uncertainties on the planet mass and radius, while Kepler-20d 
is plotted with bands illustrating the 2 a mass upper limit. Other small transiting exo- 
planets with measured masses (Kepler-lOb, CoRoT-7b, Kepler- lib cdef, Kepler-18b, 55Cnc 
e, GJ 1214b, HD 97658b, GJ 436b, Kepler-4b, HAT-P-llb) are plotted in black. The 
Solar System planets are indicated with the first letters of their names. The curves are 
illust rative consta n t-tern perature (300 K) mass-radius relations for bodies devoid of H/He 
from ISeager et al.l (120071 ). The solid lines are homogeneous-composition planets: water ice 
(blue sohd), MgSiOs perovskite (red solid), and iron (magenta solid). The non-solid lines 
are mass-radius relations for differentiated planets: 75% water ice, 22% silicate shell, and 
3% iron core (blue dashed); Ganymede-like with 45% water ice, 48.5% silicate shell, and 
6.5% iron core (blue dot-dashed); 25% water ice, 52.5% silicate shell, and 22.5% iron core 
(blue dotted); Earth-like with 67.5% silicate mantle and 32.5% iron core (red dashed); and 
Mercury-like with 30% silicate mantle and 70% iron core (red dotted). The minimal radius 
curve b ased on simulations of coUisional mantle stripping from differentiated silicate-iron 
planets ( IMarcus et al.ll2010l ) is denoted by the dashed magenta line. 
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